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MATHEMATICAL FUNCTION PACKAGE 

The following routines constitute the Mathematical Function Package 
(MFUNC), They support CALL/360-OS built-in functions and can be grouped 
as shown below. Descriptions of these routines are given on succeeding 
pages of this manual,^ in the order indicated. 

Note that, as far as possible, the general groupings are listed in 
alphabetic order according to the mnemonics of the routines. Within each 
grouping,, discussion of short floating-point precedes discussion of long 
floating-point, and discussion of real precedes discussion of complex. 

Arctangent and hyperbolic arctangent - ATAN, ATANH: 

Short Float Real Arctan (IHEATS) 

Long Float Real Arctan (IHEATL) 

Short Float Real Hyperbolic Arctan (IHEHTS) 

Long Float Real Hyperbolic Arctan (IHEHTL) 

Short Float Complex Arctan/Hyperbolic Arctan (IHEATW) 

Long Float Complex Arctan/Hyperbolic Arctan (IHEATZ) 

Error function - ERF, ERFC; 

Short Float Real Error Function (IHEEFS) 
Long Float Real Error Function (IHEEFL) 

Exponentiation - EXP: 

Short Float Real EXP (IHEEXS) 
Long Float Real EXP (IHEEXL) 
Short Float Complex EXP (IHEEXW) 
Long Float Complex EXP (IHEEXZ) 

Logarithms - LOG, LOG2, LOGIO: 

Short Float Real Log (IHELNS) 
Long Float Real Log (IHELNL) 
Short Float Complex Log (IHELNW) 
Long Float Complex Log (IHELNZ) 

Sine and cosine, hyperbolic sine and hyperbolic cosine - SIN, COS, 
SINK, COSH: 

Short Float Real Sin/Cos (IHESNS) 

Long Float Real Sin/Cos (IHESNL) 

Short Float Real Hyperbolic Sin/Cos (IHESHS) 

Long Float Real Hyperbolic Sin/Cos (IHESHL) 

Short Float Complex Sin/Cos (IHESNW) 

Long Float Complex Sin/Cos (IHESNZ) 

Square root - SQRT: 

Short Float Real SQRT (IHESQS) 
Long Float Real SQRT (IHESQL) 
Short Float Complex SQRT (IHESQW) 
Long Float Complex SQRT (IHESQZ) 

Tangent and hyperbolic tangent - TAN, TANH: 

Short Float Real Tan (IHETNS) 

Long Float Real Tan (IHETNL) 

Short Float Real Hyperbolic Tan (IHETHS) 

Long Float Real Hyperbolic Tan (IHETHL) 

Short Float Complex Tan/Hyperbolic Tan (IHETNW) 

Long Float Complex Tan/Hyperbolic Tan (IHETNZ) 



SPEED 

The average execution times given in this subsection are based on 
information in IBM Svstem/360 Instruction Timing Information (A22-6825) 
and include times for the modules called • 

ACCURACY 

In order to appreciate properly the meaning of the statistics for 
accuracy given with each module , some consideration of the limits and 
implications of these statistics is required. Because the size of 
a machine word is limited, small errors may be generated by mathematical 
routines. In an elaborate computation, slight inaccuracies can 
accumulate and become large errors. Thus, in interpreting final 
results, errors introduced during the various intermediate stages must 
be taken into account. 

The accuracy of an answer produced by a routine is influenced by two 
factors : 

1. The accuracy of the argument 

2. The performance of the routine 

Most arguments contain errors. An error in a given argument may have 
accumulated over several steps prior to the use of the routine. Even 
newly converted input data may contain slight errors. The effect of 
an argument error on the accuracy of an answer depends solely on the 
nature of the mathematical function involved and not on the particular 
coding by which that function is computed within a routine. In order 
to assist users in assessing the accumulation of errors,, a guide on 
the propagational effect of argument errors is provided for each 
function described in this subsection. Wherever possible, this guide 
is expressed as a simple formula. 

The performance statistics supplied in this document are based upon 
the assiamption that the arguments are perfect (that is, without errors 
and, therefore, having no argument error propagation effect upon 
answers). Thus, the only errors in answers are those introduced by 
the routines themselves. For each routine,, accuracy figures are given 
for the valid argument range of respective representative segments. 
In each case, the particular statistics are those most meaningful to 
the function and range under consideration. 

For example, the maximum relative error and the root-mean-square of 
the relative error of a set of answers are generally useful and 
revealing statistics. However, they are useless for the range of a 
function where its value becomes 0, since the slightest error of the 
argument value can cause an unbounded fluctuation in the relative 
magnitude of the answer. Such is the case with sin(x) for values of 
X close to pi; in this range it is more appropriate to discuss absolute 
errors. 

The statistics tabulated herein were derived from random distributions 
of 5000 arguments per segment, generated to be either uniform or 
exponential, as appropriate. It must be emphasized that each value 
quoted for the maximum error refers to a particular test using the 
method described above, and should be treated only as a guide to the 
true maximum error. 

This explains, for example, why it is possible that the maximum error 
quoted for a segment may be greater than that found from a distribution 
of different arguments over a larger range which includes the test 
range. 



HEXADECIMAL TRUNCATION ERRORS 

While the use of hexadecimal numbers in System/ 360 has led to increased 
efficiency and flexibility ^ the effect of the variable number of 
significant digits carried by the floating-point registers must be 
noted in making allowance for truncation errors. In the production 
of the CALL/360-OS PL/I library^ special care has been taken to minimize 
such errors t whenever this could be accomplished at minor cost. As 
a result g the relative errors produced by some of the library routines 
may be considerably smaller than the relative error produced in some 
instances by a single operation such as multiplication. 

Representations of finite length entail truncation errors in any number 
system. With binary normalization, the effect of truncation is roughly 
uniform. With hexadecimal normalization, however, the effect varies 
by a factor of 16 depending on the size of the mantissa; in a chain 
of computations, the worst error committed in the chain usually prevails 
at the end. 

In short- precision representation, a number has between 21 and 24 
significant binary digits. Therefore, the truncation errors range 
from 2**-24 to 2**-20 (5.96*10**-8 to 9.5*10**-7). Assuming exact 
operands,, a product or quotient is correct to the 2Uth binary digit 
of the mantissa. Hence truncation errors contributed by multiplication 
or division are no more than 2**-20. The same is true for the sum 
of two operands of the same sign. Subtraction, on the other hand, 
is the most common cause of loss of significant digits in any number 
system. For short-precision operations, therefore, a guard digit is 
provided which helps to reduce such loss. 

In long-precision representation, a number has between 53 and 56 
significant binary digits. Therefore, truncation errors range from 
2**-56 to 2**-52 C1.39*10**-17 to 2.22*10**-16) . Assuming exact 
operands,, a quotient is correct to the 56th binary digit of the 
mantissa. Therefore, truncation errors resulting from division are 
no more than 2**-52. The accuracy of a product, on the other hand, 
depends on the necessity for post-normalization. If the mantissas 
of both operands are close to 1, the truncation error of a product 
is about 2**-56. If the product of the mantissas is about 1/16, the 
truncation error is about 2**-52. On the other hand, if the mantissas 
of both operands are close to 1/16, the intermediate product has seven 
leading zeros, and post- normalization introduces four trailing zeros. 
In this case, the truncation error can be close to 2**-48 (3.5*10**-15) . 
In particular, multiplication by 1 in the long- precision form has the 
effect of erasing the last hexadecimal digit of the multiplicand. 

Normal care in numerical analysis should be exercised for addition 
and stabtraction. In particular, when two algorithms are theoretically 
equivalent, it usually pays to choose the one which avoids subtraction 
between operands of similar size. There is no guard digit for long- 
precision additions and subtractions. 

HEXADECIMAL CONSTANTS 

Many of the modules described below discriminate between algorithms 
or test for errors by comparisons involving hexadecimal constants; 
where decimal fractions are used in the descriptions, the fractions 
are only quoted as convenient approximations to the hexadecimal values 
actually employed. 



TERMINOLOGY 

Maximum and root-mean- square values for the relative and (where 
necessary) the absolute errors are given for each module. These are: 

Let f (x) = the correct value for a function 

g(x) = the result obtained from the module in question 

Then the absolute error of the result is 

ABS(f(x) - g(x)) 
and the relative error of the result is 

ABS((f (x) - g(x))/f (x)) 

Let the number of sample results obtained be N; then the root-mean- 
square of the absolute error is 

N 
SQRTI V^ (ABS(f(xi) - g (xi))**2)/N) 

I 

and the root-mean- square of the relative error is 

N 
SQRT( y^ (ABS((f(xi)-g(xl))/f(xi))**2)/N) 



ARGUMENTS 

Any restrictions on arguments are noted under two headings: 

Range : This states any range of arguments for which a module is valid. 
Arguments outside the given ranges are assumed to have been excluded 
before the module is called. 

Error and Exceptional Conditions ; These cover any conditions which 
may result from the use of a routine: 

P - Programmed conditions in the module concerned. Programmed 

tests are made where not too costly and, if an invalid argument 

is found, a branch is taken to the entry point IHEERRB of the Error 

Routine. (See "Handling of Interrupts Package" in Volume II. 

Error Routine is part of EXEP, which is a subpackage of HIP.) 

An appropriate message is printed and the ERROR condition is raised. 

I - Interrupt conditions in the module concerned. For those 
routines where SIZE and FIXEDOVERFLOW are detected by programmed 
tests or \i^ere hardware interruptions may occur, the OVERFLOW and 
UNDERFLOW conditions pass to IHEERR and are treated in the normal 
way. The machine is assumed to be enabled for all interruptions 
except significance, which is masked. 

O - Programmed conditions in modules called by the module concerned. 
These occur when invalid arguments are detected in module called. 

H - As I„ but the interrupt conditions occur in the modules called 
by the module concerned. 



TITLE: SHORT FLOAT REAL ARCTAN (IHEATS) 

Program Definition 

Piirpose and Usage 

Short Float Real Arctan is used to calculate arctan(x) or arctan(y/x) 
where x is a short floating-point real expression in radians. The 
result range is: 

arctan (x) (radians) ±pi/2 
arctan (y/x) (radians) ±pi 

Description 

Method: 

!• ATAN(y,x) 

If 

X = or 7^S(y/x) > 2**24 

the answer, SIGN(y) **pi/2g is returned except for the error 
case, X = y = 0. Otherwise 

ATAN(y,x) = ATAN(y/x)if x > or 

ATAN(y,x) = ATAN(y/x) + SIGN(y)*pi if x < 

Hence the computation is now reduced to the single argument 
case. 

2. ATAN(x) 

The general case may be reduced to the range < x < 1 since 

ATAN(-x) = - ATAN(x) and 

ATAN(l/ABS(x)) = pi/2 - ATAN(ABS(x)) 

A further reduction to the range ABS(x) < TAN (pi/12) is made 
by using 

ATAN(x) = pi/6 + ATAN((SQRT(3)*x-l)/(x ♦ SQRTO))) 

Care is taken to avoid the loss of significant digits in 
computing SQRTO) ♦x - 1. 

For the basic range ABS(x) < TAN (pi/12) use an approximation 
formula of the form 

ATAN(x)/x = a ♦ b*x**2 + c/ (d + x**2) 

with relative error less than 2**-27.1. 

Effect of an Argument Error : 

Let 

t = X or y/x 

Then the absolute error of the answer approaches the absolute error 
in t divided by (1 + t**2) . Hence, for small values of t, the two 



errors are approximately the same; however, as t becomes larger, the 
effect of the argument error on the answer error diminishes. 

Acciaracy : 



j Arguments | 

I I 

, ^. 

I Range | Distribution | 

L ^ 



Relative Error 
♦106 



RMS 



I 



Maximum 



IHEATSl 



Full 
Range 



IHEATS2 



Full 
Range 



Tangents of 

numbers 

uniformly 

distributed 

between 

-pi/2 and 

pi/2 



0.443 



0.958 



y = sin and 
X = cos of 
numbers 
uniformly 
distributed 
between 
-pi/2 and 
pi/2 



0.449 



1.42 



Implementation ; 

• Module size: 408 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the appropriate entry point 
in the tables: 

k = TANDCIB) 



IBM System/360 Model Number 
I 40 I 50 I 65 I 



I ABS(x) 

L 



30 



75 



IHEATSl 



1 ^k 1 

1 +- 

1 k < 1 
1 ABS(x) 1 

1 <1 1 
1 +- 

1 1^ 1 
1 ABS(x) 1 

1 <l/k 1 
1 1_ 

1 ^1/k 1 


3162 1 858 1 

1 1 
U3a5 1 1136 1 

1 1 
^ +_ 

1 1 
5051 1 1301 1 

1 1 

+ ^.- 

3868 1 1023 1 


279 1 78. a 1 

1 1 
3«»7 1 97.8 1 

1 1 

+ ^.- 

1 1 

381 1 108 1 
i 1 

313 1 88.7 1 


1 

48. It 1 
+ 

58.0 1 


64.1 1 


54.5 1 



IHEATS2 

^ <k I 4193 I 1138 I 363 | 106 I 67.3 

k< I I I I I 

ABS(x) I 5376 | 1416 | 431 | 126 | 77.0 

<1 I I I I I 
1 ^ ^ 1 1 

1^ I I I I I 

ABS(x) I 6082 I 1581 | 465 | 136 | 83.1 

<l/k ^ 1 ^ ' ^ I I I ^ 

^1/k j 4899 I 1303 | 496 | 117 | 73.4 



Errors Detected 

Error and Exceptional Conditions ; 

P: IHEATS22 x = y = (206) 
Local Variables 
None 

Program Interface 
Entry Points 

IHEATSl - Entry for ATAN(x) 



P7 
PLIST 



= A (PLIST) 
= A(x) 

A (target) 



IHEATS2 = Entry for ATAN(y,x) 



P7 
PLIST 



= A (PLIST) 
= A(y) 

A(x) 

A (target) 



Exit Conditions 

Normal. Returns to caller via the link register. 
Abnormal. Branches to EXEP via the entry point IHEERRB. 



Routines Called 

EXEP Execution Error Package 

Global Variables 

None 

Coinroents 

IHEATSl is called by compiled code, IHEATS2 is called by compiled 
code, IHELNW, and IHEATW, 



TITLE: LONG FLOAT REAL ARCTAN (IHEATL) 

Program Defini1:ion 

Purpose and Usage 

Long Float Real Arctan is used to calculate arctan(x) or arctan(y/x), 
where x is a long floating-point real expression in radians • The 
result range is: 

arctan (x) (radians) ±pi/2 
arctan (y/x) (radians) ±pi 

Description 

Method: 

1. ATAN(y,x) 

If X = or ABS(y/x) > 2**56, the answer SIGN(y)*pi/2 is returned 
except for the error case x = y = 0, Otherwise 

ATAN(y,x) = ATAN(y/x) if x > or 

ATAN(y,x) = ATAN(y/x) + SIGN(y)*pi if x < 

Hence the computation is now reduced to the single argument 
case. 

2. ATAN(x) 

The general case may be reduced to the range < x < 1 since 

ATAN(-x) = - ATAN(x) and 

(ATAN(l/ABS(x)) = pi/2 - ATAN(ABS(x)) 

A further reduction to the range ABS(x) < TAN(pi/12) is made 
by using 

ATAN(x) = pi/6 + ATAN((SQRT(3)*x-l)/(x + SQRTO))) 

Care is taken to avoid the loss of significant digits in 
computing 

SQRT(3)*x - 1 

For the basic range ABS(x) < TAN (pi/12), use a continued fraction 
of the form 

ATAN(x)/x = 1 + al*x*x/(bl + x*x + a2/(b2+x*x + a3/(b3 + x*x 
+ a4/(b4 ♦ x*x)))) 

with relative error less than 2**(-57,9). 

Effect of an Argument Error : 

Let 

t = X or y/x 

Then the absolute error of the answer approaches the absolute error 
in t divided by (1 + t**2). Hence, for small values of t, the two 
errors are approximately the same; however, as t becomes larger, the 
effect of the argument error on the answer error diminishes. 



Accuracy : 

IHEATLl 



Arguments 



1 
I 



Relative Error 
♦10^5 



Range | Distribution! RMS | Maximum 
1 ^ ^ 

-l<x<l I Uniform | 0.0438 | 0,207 



Implementation ; 

• Module size: 544 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the appropriate entry point 
in the table: 



r 

1 

1 _ _ 


1 




IBM Systero/360 Model Number 




1 Entry 

1 Point 
, 

1 IHEATLl 
1 IHEATL2 


1 

1 

-+- 

1 

1 


30 

20472 
23523 


1 40 1 50 1 65 1 

1 1 1 1 
_^. ^ ^ ^._ 

1 4389 1 826 1 181 | 
1 5042 1 967 1 217 | 


75 

100 
123 



Errors Detected 

Error and Exceptional Conditions ; 

P: IHEATL2: x = y = (206) 
Local Variables 
None 

Program Interface 
Entry Points 

IHEATLl - Entry for ATAN(x) 

P7 = A(PLIST) 
PLIST = A(x) 

A (target) 

IHEATL2 - Entry for ATAN(y,x) 

P7 = A (PLIST) 
PLIST = A(y) 

A(x) 

A (target) 

Exit Conditions 

Noinnal, Returns to caller via the link register • 
Abnormal. Branches to EXEP via the entry point IHEERRB. 



10 



Routines Called 

EXEP Execution Error Package 

Global Variables 

None 

Comments 

IHEATLl is called by compiled code. IHEATL2 is called by compiled 
code, IHELNZ, and IHEATZ, 



11 



TITLE; SHORT FLOAT REAL HYPERBOLIC ARCTAN (IHEHTS) 

Program Definition 

Purpose and Usage 

Short Float Real Hyperbolic Arctan is used to calculate hyperbolic 
arctanCx)^ where x is a short floating-point real expression in radians. 

Description 

Method ; 

1. ABS(x) < 0. 2 

Use a rational approximation of the form: 
ATANH(x) = X ♦ x**3/(a + b*x**2) 

2. 0^2 < ABS(x) < 1 

ATANH(x) = -SIGN(x)*0.5*LOG((0.5 - ABS{x/2) )/(0.5 + ABS(x/2))) 

Effect of ^ Argument Error ; 

The absolute error caused in the result is approximately equal to the 
absolute error in the argument divided by (1 - x**2) . Thus as x 
approaches +1 or -1, relative error increases rapidly. Near x = 0, 
the relative error in the result is of the order of that in the 
argument. 

Accuracy ; 



Arguments 


1 
1 


Relative 
*10e 


Error 




J 










T 








Range | Distribution j 


RMS 


1 


Maximum 


^ 


-+- 




-+ 




-0.8<x 1 Uniform 


1 


o.«»«»o 


1 


1.32 


<0.8 1 


1 




1 




+ 


-+- 








-0.9<x 1 Uniform 


1 


0.389 


1 


I.IU 


<0.9 1 


1 




1 





Implementation ; 

• Module size; 192 bytes 

• Execution times; 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the table; 
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r - - , 

I I IBM Syst:em/360 Model Number 

I 

I ABS(x) I 30 I UO I 50 I 65 I 75 

, 

|<0.2 I 2520 I 667 | 208 | 52.0 | 31.3 

|0.2< I I I I I 
JABSCx) I 7091 I 1829 | 606 | 163 | 94.8 

I<1 I I I I I 

t 

Errors Detected 

Error and Exceptional Conditions ; 

P: ABS(x) > 1 (208) 
Local Variables 
None 

Program Interface 
Entry Points 

IHEAHSO 

P7 = A(PLIST) 
PLIST = A(x) 

A (target) 

Exit Conditions 

Normal. Returns to caller via the link register. 
Abnormal. Branches to EXEP via the entry point IHEERRB. 

Routines Called 

IHELNS Short Float Real Log 
EXEP Execution Error Package 

Global Variables 

None 

comments 

Called by compiled code and IHEATW. 
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TITLE: LONG FLOAT REAL HYPERBOLIC ARCTAN (IHEHTL) 

Program Definition 

Purpose and Usage 

Long Float Real Hyperbolic Arctan is used to calculate hyperbolic 
arctan(x), where x is a long floating-point real expression in radians, 

Description 

Method: 

1. ABS(x) < 0.25 

Use a Chebyshev polynomial of degree 8 in x**2 to compute 
ATANH(x)/x. 

2. 0,25 < ABS(x) < 1 

ATANH(x) = -SIGN(x)*0.5*LOG((0.5 - ABS(x/2) )/(0.5 + ABS(x/2))) 

Effect of an Argument Error ; 

The absolute error caused in the result is approximately equal to the 
absolute error in the argument divided by (1 - x**2) . Thus as x 
approaches +1 or -1, relative error increases rapidly. Near x = 0, 
the relative error in the result is of the order of that in the 
argument. 

Accuracy : 



I Arguments | Relative Error 

I I *10^5 

, 1 

I Range | Distribution! RMS | Maximum 

|ABS(x)<| Uniform | 0.0650 | 0.223 
10.25 I ^ 1^ I 

|ABS(x)<| Uniform j 0.133 | 0.397 
10.95 I I I 

L J 

Implementation ; 

• Module size: 262 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the table: 
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I I IBM System/360 Model Number 

I 

I ABS(x) I 30 I 40 I 50 I 65 I 75 



|<0,25 I 12252 I 3037 | 562 | 121 I 68 

|0.25< I I I i I 
|ABS(x) I 20448 | 4900 | 1040 1 242 | 137 

I<1 I I I I I 

L , 

Errors Detected 

Error and Exceptional Conditions ; 

P: ABS(x) > 1 (208) 
Local Variables 
None 

Program Interface 
Entry Points 

IHEAHLO 

P7 = A(PLIST) 
PLIST = A(x) 

A (target) 

Exit Conditions 

Normal, Returns to caller via the link register. 
Abnormal. Branches to EXEP via the entry point IHEERRB. 

Routines Called 

IHELNL Long Float Real Log 
EXEP Execution Error Package 

Global Variables 

None 

Comments 

Called by compiled code and IHEATZ. 
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TITLE: SHORT FLOAT COMPLEX ARCTAN/HYPERBOLIC ARCTAN (IHEATW) 

Program Definition 

Purpose and Usage 

Short Float Complex Arctan/Hyperbolic Arctan is used to calculate 
arctan(z) or hyperbolic arctan (z), where z is a short floating-point 
complex expression in radians. 

Description 

Method: 

Let 

z = X + yl 
Then 

REAL(ATANH(z)) = (ATANH (2*x/(l ♦ x*x ♦ y*y)))/2 

IMAGCATANHCz)) = (ATAN(2*y, (1 - x*x - y*y)))/2 and 

ATAN(z) = -(ATANH(zI))I 

Effect of an Argument Error: 

The absolute error in the result is approximately equal to the absolute 
error in the argument divided by (1 ♦ z**2) in the case of IHEATTO, 
or by (1 - z**2) in the case of IHEAHTO. Thus the effect may be 
considerable in the vicinity of z = IIKIHEATTO) or 
z = ±1 (IHEAHTO). 

Accuracy : 

r T 

I Arguments | Relative Error | 

I I *10« I 

I 1 , 

I Range I Distribution | RMS | Maximum | 
IHEATTO 

fruii I 11 1 

I Range I Exponential | 0.216 | 2.88 | 

L J 

IHEAHTO 

[fuii I II 1 

I Range I Exponential | 0.208 | 1.18 | 
L J 

Implementation ; 

• Module size: 304 bytes 

• Execution times: 
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Approximate execution times in microseconds for the System/360 
models given below are obtained from the appropriate entry point 
in the table ^ where 

a = ABS(2*u/(l+x*x+y*y)) 
u = y for IHEATTO 
= x for IHEAHTO 

I I IBM System/ 360 Model Number | 
, , 

I I 30 I 40 I 50 I 65 I 75 I 

L J 

IHEATTO 

^ ^ 

|a<0-2 I 12235 | 3306 | 1033 | 279 | 173 | 

|0.2<a I 16056 | «»45a | 1«*08 | 462 | 276 | 
I<1 I I I I I I 

IHEAHTO 

|a<0,2 I 12100 I 3267 | 1017 | 275 | 171 | 

|0.2<a I 15921 | 4415 | 1392 | 458 | 274 | 
I<1 I I I I I I 

L J 

Errors Detected 

Error and Exceptional Conditions ; 

P: IHEATTO z = ±11 (211) 
IHEAHTO z = ±1 (211) 

Local Variables 

None 

Program Interface 

Entry Points 

IHEATTO - Entry for ATAN(z) 

P7 = A(PLIST) 
PLIST = A(z) 

A (target) 

IHEAHTO - Entry for ATANH(z) 

Linkage as for IHEATTO 

Exit Conditions 

Normal. Returns to caller via link register. 

Abnormal. Branches to EXEP via the entry point IHEERRB. 



17 



Routines Called 

IHEATS Short Float Real Arctan 

IHEHTS Short Float Real Hyperbolic Arctan 

EXEP Execution Error Package 

Global Variables 

None 

Comments 

Called by compiled code^ 
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TITLE: LONG FLOAT COMPLEX ARCT AN/HYPERBOLIC ARCTAN (IHEATZ) 

Program Definition 

Purpose and Usage 

Long Float Complex Arctan/Hyperbolic Arctan is used to calculate arctan(z) 
or hyperbolic arctan (z),, where z is a long floating-point complex 
expression in radians • 

Description 

Method ; 

Let 

z = X ♦ yl 
Then 

REAL(ATANH(z)) = (ATANH(2*x/(l + x*x ♦ y*y)))/2 

IMAG(ATANH(z)) = (ATAN(2*y,, (1 - x*x - y*y)))/2 and 

ATAN(z) - -(ATANH(zI))I 

Effect of an Argument Error : 

The absolute error in the result is approximately equal to the absolute 
error in the argument divided by (1 + z**2) in the case of IHEATMO, 
or by (1 - z**2) in the case of IHEAHMO, Thus the effect may be 
considerable in the vicinity of z = ilKIHEATMO) or 
z = IKIHEAHMO). 

Accuracy : 

I Arguments | Relative Error | 
I I *10^5 I 

, + , 

I Range | Distribution) RMS | Maximum | 
L J 

IHEAIMO 

^ ^ 

I Full (Exponential | 0.141 | 7.93 | 
I Range | I | I 

L J 

IHEAHMO 

I Full I Exponential | 0.0826 | 1.20 | 
I Range | I | I 

L J 

Implementation ; 

• Module size: 296 bytes 

• Execution times: 
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Approximate execution times in microseconds for the System/360 
models given below are obtained from the appropriate entry point 
in the table, where 

a = ABS( 2*11/(1 ♦ x*x + y*y)) 

u = y for IHEATMO 
= X for IHEAHMO 

^ , ^ 

I I IBM System/ 360 Model Number | 
I , 

I I 30 I UO I 50 I 65 I 75 I 

L J 

IHEATMO 

|a<0.25 I i\3iM7 \ 9977 | 2006 | 455 | 260 1 

|0,25<a j 51673 | 11840 | 2406 | 576 | 329 | 
I<1 I I I I II 

L J 

IHEAHMO 

r : T 

|a<0,25 I 43293 | 9923 | 1987 | 450 | 258 | 

|0.25<a I 51489 | 11786 | 2466 | 571 | 327 | 
I<1 I I I I II 

Errors Detected 

Error and Exceptional Conditions ; 

P: IHEATMO z = ±11 (211) 
IHEAHMO z = ±1 (211) 

Local Variables 

None 

Program Interface 

Entry Points 

IHEATMO - Entry for ATAN(z) 

P7 = A(PLIST) 
PLIST = A(z) 

A (target) 

IHEAHMO - Entry for ATANH(z) 

Linkage as shown for IHEATMO 

Exit Conditions 

Normal, Returns to caller via the link register. 
Abnormal, Branches to EXEP via the entry point IHEERRB, 
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Routines Called 

IHEATL Long Float Real Arctan 

IHEHTL Long Float Real Hyperbolic Arctan 

EXEP Execution Error Package 

Global Variables 

None 

Comments 

Called by compiled code. 



21 



TITLE: SHORT FLOAT REAL ERROR FUNCTION (IHEEFS) 

Program Definitdon 

Purpose and Usage 

Short Float Real Error Function is used to calculate the error function 
of X or the complement of this function,, \^ere x is a short floating- 
point real expression. 

Description 

Method ; 

1. < X < 1.317 

Compute ERF(x)/x by using a Chebyshev interpolation polynomial 
of degree 6 in x**2, with relative error less than 2**-24. 

ERFCCx) = 1 - ERF(x) 

(ERFC(x) > 1/16 in this range.) 

2. 1.317 < X < k, where k = 2.0«*000092 

compute ERFCCx) by using a Chebyshev interpolation polynomial 
of degree 7 in (x-k) , with absolute error less than 1.3*2**-30. 

ERFCx) = 1 - ERFCCx) 

CERFCCx) > 1/256 in this range.) 

3. k < X < 13.306 

ERFCCx) ♦xEXPCx**2) is computed by using a Chebyshev interpolation 
polynomial of degree 6 in x**-2, with relative error less than 
1^2*2**-23: 

If X < 3.9192, ERFCx) = 1 - ERFCCx) 

If X > 3.9192, ERFCx) = 1 

4. X > 13.306 

Results 1 and are returned for ERFCx) and ERFCCx) respectively. 

5. X < 

ERFCx) = -ERFC-x) 

Effect of an Argument Error ; 

The absolute error caused in the result is approximately equal to the 
absolute error in the argument multiplied by EXPC-x**2). 

ERFCx) ; As the magnitude of the argument increases from 1, the effect 
of an argument error diminishes rapidly. For small x, the relative 
error of the result is of the order of the relative error of the 
argument. 
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Accuracy ; 



.——-■'——— 


Arguments 


1 
1 


Relative 
♦10« 


Error 


Range 


- T 
1 Distribution | 


RMS 1 


Maximum 


ABS(x) 


1 Unifoinn 




0.139 1 


0.934 


<1,3 


-X —.-.——.-.-.——. 


-.X. 


_____!___ 




1.3< 


-^ 


T 


+ 




ABS(x) 


1 Uniform 




0-0372 1 


0.263 


<2 










2< 


-+ 


-+- 


1 




ABS (x) 


1 Uniform 




0.03U7 1 


0.0605 


<3.9 











Implementation ; 

• Module size; 376 bytes 

• Execution times; 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the appropriate entry point 
in the table; 



IBM System/360 Model Number 
I 40 I 50 I 65 I 



|ABS(x) 

L 



30 



75 



IHEEFSO 



|<1.32 

|1.32< 
|ABS(x) 

|<2.04 

(2,04< 

|ABS(x) 

|<3.92 

|>3,92 
L 



I 4354 I 1191 I 392 

I I i 

I 4613 I 1266 I 418 

I I I 

I I I 

I 10013 I 2843 I 868 

I II 

I 1530 I 473 I 183 



I 100 I 57,2 

I I 

I 110 I 62.0 

I I 

I I 

I 228 I 140 

I I 

I 50,7 I 32.0 



Errors Detected 

None 

Local Variables 

None 
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Program Interface 

Entry Points 

IHEEFSO ~ Entry for ERF(x) 

P7 = A(PLIST) 
PLIST = A(x) 

A(Target) 

Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

IHEEXS Short Float Real EXP 

Global Variables 

None 

Comments 

Called by compiled code. 
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TITLE: LONG FLOAT REAL ERROR FUNCTION (IHEEFL) 

Program Definition 

Purpose and Usage 

Long Float Real Error Function is used to calculate the error function 
of X or the complement of this function, where x is a long floating- 
point real expression. 

Description 

Method: 

1. < X < 1 

Compute ERF(x>/x by using a Chebyshev interpolation polynomial 
of degree 11 in x**2, with relative error less than 1.07*2**-57. 

ERFC(x) = 1 - ERF(x) 

(ERFC(x) > 1/16 in this range) 

2. 1 < X < 2.0U000092 

Compute ERFC(x) by using a Chebyshev interpolation polynomial of 
degree 18 in (x - !• 999999), with absolute error less than 
1.5*2**-61. 

ERF(x) = 1 - ERFC(x) 

(ERFC(x) > 1/256 in this range) 

3. 2,04000092 < x < 13,306 

ERFC(x) is computed by using a Chebyshev interpolation polynomial 
of degree 20 in x**2 for ERFC(x)*xEXP(x**2) , with relative error 
ranging from 2**-53 at 2,04000092 to 2**-51 at 13.306: 

If X < 6.092, ERF(x) = 1 - ERFC(x) 

If X > 6,092, ERF(x) = 1 

4. X > 13,306 

Results 1 and are returned for ERF(x) and ERFC(x), respectively. 

5. X < 

ERF(x) = -ERF(-x) 

Effect of an Argument Error : 

The absolute error caused in the result is approximately equal to the 
absolute error in the argument multiplied by EXP(-x**2), 

ERF(x): As the magnitude of the argument increases from 1, the effect 
of an argument error diminishes rapidly. For small x, the relative 
error of the result is of the order of the relative error of the 
argument. 
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Accuracy ; 



Arguments 



j Range | Distribution 

L 



|ABS(x)<|Unifonn 
11-317 I 

I 



|2-04< lUniforin 
|ABS(x) I 

|<6.092 I 

L 



Relative Error 
♦10^5 



RMS 
0.0280 



I 



0.0107 



|1.317< I Uniform 
|ABS(x) I 
|<2.04 I 

0.00803 



Maximum 
0.202 



0.0291 



0.0170 



Implementation ; 

• Module size: 7«*U bytes 

• Execution tiroes: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the appropriate entry point 
in the table: 



IBM System/ 360 Model Number 



|ABS(x) 
I 



I 



30 



I 



i\0 



I 



50 



I 



65 



I 



75 



|<1.00 I 16567 I iH54 | 805 

|1.00< I I I 

|ABS(x) I 24413 I 6095 | 1175 

|2.04< I I I 

|ABS(x) I 45574 I 1105 | 2142 

l<6^09 I 1^ ^ 1^^ 

|>6.09 I 2493 I 707 \ 222 

t 



I 180 I 103 

I I 

I 263 I 147 

I ^ I 

I I 

I 477 I 269 

I ^ 1^ 

I 58.0 I 36.3 



Errors Detected 

None 

Local Variables 

None 

Program Interface 

Entry Points 

IHEEFLO - Entry for ERF(x) 

P7 = A(PLIST) 
PLIST = A(x) 

A (target) 
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Exit Conditions 

Normal. Returns to caller via the link register • 

Routines Called 

IHEEXL Long Float Real EXP 
Global Variables 
None 
Comments 
Called by compiled code. 
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TITLE: SHORT FLOAT REAL EXP CIHEEXS) 

Program Definition 

Purpose and Usage 

Short Float Real EXP is used to compute e to the power x,, where x is a 
short floating-point real expression. 

Description 

Method ; 

If x < -180. 2, a zero result is returned immediately. Otherwise, 
EXP(x) is calculated as 2**(x*LOG2(e) ) . 

The calculation is performed as follows: 

x*LOG2(e) = f ♦ N 

where N = 4h + g, h is an integer such that g = 0, 1, 2, or 3, and 
< f < 1.0. 

Then, by subtracting 0.5,, this is reduced to the range -0.5 < f < 0.5, 
Next, 2**f is calculated as (a ♦ b/(c + x*x) + x)/(a + b/(c -i^ x*x) - x) . 

This is multiplied by 2**0.5 and then shifted in the appropriate 
direction to give the effect of multiplication by 2**g. Finally, 
the exponent of the result is obtained from h. 

Effect of an Argument Error: 

The relative error caused in the result is approximately equal to the 
absolute error in the argument, that is,, to the argument relative error 
multiplied by x.. Thus, for large values of x, even the roundoff error 
of the argument causes a substantial relative error in the answer. 

Accuracy : 



Arguments 



I 
I 



Relative Error 
*10« 



Range | Distribution! RMS | Maximum 



-l<x<l I Uniform 

Full I Uniform 
Range | 



-+ +- 

I 0.132 I 0.490 

I 1.29 I 2.61 
I I 



Implementation : 

• Module size: 232 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the table: 
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r 


1 




IBM Systein/360 Model Number 






1 


30 


t 40 1 50 1 65 1 


75 


L 


1 


38U7 


1 1172 1 356 1 90.0 | 


58.0 



Errors Detected 

Error and Exceptional Conditions ; 

I: OVERFLOW if x > 17«».673 (201) 
Local Variables 
None 

Program Interface 
Entry Points 

IHEEXSO 

P7 = ACPLIST) 
PLIST = A(x) 

A (target) 

Exit Conditions 

Normal • Returns to caller via the link register. 

Routines Called 

None 

Global Variables 

None 

Comments 

Called by compiled code, IHEXXS, IHESHS, IHETHS, IHEEFS, IHEEXW, and 
IHESNW. 
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TITLE: LONG FLOAT REAL EXP (IHEEXL) 

Program Definit:ion 

Purpose and Usage 

Long Float Real EXP is used to compute e to the power x,, where x is a 
long floating-point real expression. 

Description 

Method ; 

If X < -130.2183^ return zero result. Otherwise, let 

y = x/L0G(2) 

= 4*a - b - c/16 - d 

where a, h, and c are integers, 0<b<3, 0<c< 15, and < d < 1/16, 

Then 

EXP(x) = 2**y 

= 16**a*2**(-b)*2**(-c/16)*2**(-d) • 

Compute 2**(-d) by using the Chebyshev interpolation polynomial of 
degree 6 over the range < d < 1/16, with maximum relative error 
2**(-57). If c > 0, multiply 2** (-d) by 2** (-c/16). The constants 
2**(-c/16), 1 < c < 15„ are included in the subroutine. If b > 0, 
halve the result b times. Finally, multiply by 16**a by adding a to 
the characteristic of the result. 

Effect of an Argument Error ; 

The relative error caused in the result is approximately equal to the 
absolute error in the argument, that is, to the argument relative error 
multiplied by x. Thus, for large values of x, even the roundoff error 
of the argument causes a substantial relative error in the answer. 

Accuracy ; 



1 
1 


Arguments 


1 

1 

1 


Relative Error | 
♦10" 1 


1 

1 Range 
1 «— — -.— 


"" T 

1 Distribution | 


RMS 1 


1 

Maximum | 


|-1<X<1 
1 — 


1 Uniform 


T 
1 


0.0674 1 


0.216 1 


IFull 
I Range 


1 Uniform 

1 


T 

1 
1 


0.867 1 
1 


2.30 1 



Implementation ; 

• Module size; 448 bytes 

• Execution times; 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the table; 
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r 


1 


IBM Systein/360 Model Number 






1 30 


1 i>0 1 50 1 65 i 


75 


L 


1 12131 


1 2901 1 616 1 343 | 


194 



Errors Detected 

Error and Exceptional Conditions ; 

I: OVERFLOW if x > 17U.673 (201) 
Local Variables 
None 

Program Interface 
Entry Points 

IHEEXLO 

P7 = A(PLIST) 
PLIST = A(x) 

= A (target) 

Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

None 

Global Variables 

None 

Comments 

Called by compiled code, IHEXXL, IHESHL, IHETHL, IHEEFL, IHEEXZ, and 
IHESNZ . 
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TITLE: SHORT FLOAT COMPLEX EXP (IHEEXW) 

Program Definition 

Purpose and Usage 

Short Float Complex EXP is used to calculate e to the power z, where z 
is a short floating-point complex expression. 

Description 

Method ; 

Let 

z = X + yl 
Then 

REAL (EXP (z)) = EXP(x)*COS(y) and 

IMAG(EXP(z)) = EXP(x)*SIN(y) 
Effect of an Argument Error ; 
Let 

EXP(x ♦ yl) = s*EXP(kI) 

Then k = y, and the relative error in s is approximately equal to the 
absolute error in x. 

Accuracy ; 



r 












1 


1 


Arguments 


1 


Relative 


i Error 




1 

1 . • ^. 




1 




♦10« 






1 




— t 










1 Range 


1 Distribution | 


RMS 


1 


Maximum 




1 





~+- 




•~+ 




— 1 


|ABS(x) 


1 Uniform 




1.32 




3.14 




|<170 














|ABS(y) 














|<pi/2 

1 —-.—-.—.- 




^.j- 




J- — 






1 


■+ ----- 


T 




T 






|ABS(x) 


1 Uniform 




1.31 




3.34 




|<170 














|pi/2< 














|ABS(y) 














|<20 

L 












— J 



Implementation ; 

• Module size: 136 bytes 

• Execution tiroes; 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the table: 
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r 


1 


IBM System/360 Model Number 






1 30 


1 40 1 50 1 65 1 


75 




1 14240 


1 4058 1 1244 1 326 | 


200 



Errors Detected 

Error and Exceptional Conditions ; 

O: ABS(y) > 2**18*pi (212) in Short Float Real Sin/Cos (IHESNS) 
H: OVERFLOW (300) in Short Float Real EXP (IHEEXS) 

Local Variables 

None 

Program Interface 

Entry Points 

IHEEXTO 

P7 = A(PLIST) 
PLIST = A(z) 

A (target) 

Exit Conditions 

Normal • Returns to caller via the link register. 

Routines Called 



IHESNS Short Float Real Sin/Cos 
IHEEXS Short Float Real EXP 

Global Variables 

None 

Comments 

Called by compiled code and IHEXXW. 
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TITLE: LONG FLOAT COMPLEX EXP (IHEEXZ) 

Program Definition 

Purpose and Usage 

Long Float Complex EXP is used to calculate e to the power z, where z 
is a long floating-point complex expression • 

Description 

Method: 

Let 

z = X + yl 
Then 

REAL (EXP (z)) = EXP(x)*COS(y) and 

IMAG{EXP(z)) = EXP(x)*SIN(y) 
Effect of an Argument Error : 
Let 

EXP(x + yl) = s*EXP(kI) 

Then k = y, and the relative error in s is approximately equal to the 
absolute error in x. 

Accuracy : 



r 
1 
1 

1 — «—- — .- 


Arguments | 


Relative Error 
♦10" 


\ 

1 Range 

1 

jABSCx) 

I<1 
|ABS(y) 

|<pi/2 

1 


— ~ "t 

1 Distribution | 

1 Uniform | 
-I 1- 


RMS 1 Maximum 
^ 

0.136 1 0.«*78 
^ 


|ABS(x) 

|<20 

jABSCy) 

|<20 

1 


1 Uniform | 


1.28 1 2.29 



Implementation : 

• Module size: 136 bytes 

• Execution times: 

Approximate execution times in microseconds for System/360 models given 
below are obtained from the table: 
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r 


1 


IBM System/360 Model Number 






1 


30 1 aO 1 50 1 65 1 


75 


L 


1 


U2838 1 10560 | 217i» | 505 | 


287 



Errors Detected 

Error and Exceptional Conditions ; 

O: ABS(y) > 2**50*pi (203) in Long Float Real Sin/Cos (IHESNL) 
H: OVERFLOW (300) in Long Float Real EXP (IHEEXL) 

Local Variables 

None 

Program Interface 

Entry Points 

IHEEXMO 

P7 = A(PLIST) 
PLIST = A(z) 

A (target) 

Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

IHESNL Long Float Real Sin/Cos 
IHEEXL Long Float Real EXP 

Global Variables 

None 

Comments 

Called by compiled code and IHEEXL. 
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TITLE: SHORT FLOAT REAL LOG (IHELNS) 

Program Definition 

Purpose and Usage 

Short Float Real Log is used to compute: 

Log X to the base e LOG(x) 
Log X to the base 2 L0G2(x) 
Log X to the base 10 LOGlO(x) 

where x is a short floating-point real expression. 

Description 

Method ; 

Let 

X = m*16**p 

where 1/16 < m < 1 and p is an integer. 

Two constants^ a (= base point) and b (= -LOG2(a))^ are defined as 
follows : 



1/16 


< m < 1/8: 


a = 1/16 


b = 4 


1/8 


< m < 1/2: 


a = 1/4 


b = 2 


1/2 


< m < 1: 


a = 1 


b = 



Let 

y = (m - a)/(m ♦ a) 
Then 

m = a*(l+y)/(l-y) and ABS(y) < 1/3 
Now 

X = 2**(4*p-b)*(l+y)/(l-y) 
Therefore 

LOG(x) = («»*p-b)*LOG(2) + LOG( (l+y)/(l-y) ) 

LOG( (l+y)/(l-y)) is computed by using the Chebyshev interpolation 
polynominal of degree 4 in y**2 for the range <y**2<l/9, with 
maximum relative error 2**-27.8. LOG2 (x) or LOGIO (x> is calculated 
by multiplying LOG(x) by LOG2(e) or LOGIO (e) respectively. 

Effect of an Argument Error : 

The absolute error caused in the result is approximately equal to the 
relative error in the argument. Thus if the argument is close to 1, 
even the roundoff error of the argument causes a substantial relative 
error in the answer j, since the function value there is very small. 
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Accuracy ; 



Arguments 



I 

I i 

I 1, 

I Range | Distribution | 

L 



Relative Error 
*10« 



RMS 



Maximum 



IHELNSO 



I Exclud- I I I 

ling I I I 

|0.5<x I Exponential] 0.032 | 

l<2.0 I I I 

L — 



0,577 



IHEL2S0 



j Exclud- I 
ling I 
|0,5<x I 
|<2.0 I 

L 



I 



Exponential] 0.3U2 ] 0,75t* 



I 



I 



IHELGSO 



] Exclud- I ] 

ling 1 I 

] • 5<x I Exponential J 



I 
I 

0.170 ] 0.963 



|<2.0 

L_ 


1 


1 


1 






r " 

1 
1 

1 


Arguments 


1 

1 


Absolute Error 
*10« 


1 

] Range 


T 

1 Distribution ] 


RMS 1 


Maximum 


IHELNSO 


|0.5<x 
]<2.0 
L— -> — — 


1 Uniform 
1 


1 
1 


0.0960 1 

1 


0.39^ 


IHEL2S0 


10.5<x 
]<2.0 


] Uniform 

1 


1 
1 


0.177 1 

1 


0.8U2 


IHELGSO 


10.5<x 

|<2.0 

I 


1 Uniform 

1 


1 
1 


0.0526 I 

1 


0.164 
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Implementation ; 

• Module size: 256 bytes 

• Execution tiroes: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the appropriate entry point 
in the table: 

^ ^ 

I I IBM System/ 360 Model Number 

, 



(Entry | 30 | 
1 Point 1 1 


40 1 50 


1 65 1 75 
1 1 


1 — 

IIHELNSO 1 4669 | 

, 1 ^. 

IIHEL2S0 1 5041 | 
1 ^ 1. 

IIHELGSO 1 5054 | 

L 


1238 1 385 
+ 

1342 1 417 

1366 1 417 


1 173 I 95.7 
.«+ + 

1 180 1 101.3 
..^ ^ 

1 180 1 101.3 



Errors Detected 

Error and Exceptional Conditions : 

P: X < (202) 
Local Variables 
None 

Program Interface 
Entry Points 

IHELNSO - Entry for LOG(x) 



P7 


= 


A (PLIST) 


PLIST 


~ 


A(x) 

A (Target) 



IHELGSO - Entry for LOGlO(x) 

Linkage as for IHELNSO 
IHEL2S0 - Entry for L0G2(x) 

Linkage as for IHELNSO 

Exit Conditions 

Normal. Returns to caller via the link register. 
Abnormal. Branches to EXEP via the entry point IHEERRB. 

Routines Called 

EXEP Execution Error Package 

Global Variables: None 

Comments 

Called by compiled code,, IHEXXS, IHEXXW, IHEHTS, and IHELNW, 
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TITLE: LONG FLOAT REAL LOG (IHELNL) 

Program Def inii:ion 

Purpose and Usage 

Long Float Real Log is used to compute: 

Log X to the base e LOG(x) 
Log X to the base 2 LOG2(x) 
Log X to the base 10 LOGlO(x) 

where x is a long floating-point real expression • 

Description 

Method ; 

Let 

X = 16**p*2**(-g)*ra 

where p is the exponent, q is an integer such that < q < 3, 
and 1/2 < m < 1. 

Two constants r a (= -base point) and b (= LOG2(a)), are defined 
as follows: 

1/2 < m < 1/SQRT(2): a = 1/2, b = 1 

1/SQRT(2) <m<l: a=l, b=0 
Let 

y = (m - a)/(m ••• a) 
Then 

m = a*(l + y)/(l - y) 
and 

ABS(y) < 0.1716 
Now 

X = 2**(4*p - q - b)*(l + y)/(l - y) 
Therefore 

LOG(x) = (U*p - q - b)*L0G(2) ♦ LOG((l + y)/(l - y) ) 

L0G((1 + y)/(l - y)) is computed by using the Chebyshev interpolation 
polynominal of degree 7 in y**2 for the range < y**2 < 0.029U4, with 
maximum relative error 2** (-59. 6). 

L0G2(x) or LOGlO(x) is calculated by multiplying the result by L0G2(e) 
or LOGlO(e) respectively, 

Effect of an Argument Error: 

The absolute error caused in the result is approximately equal to the 
relative error in the argument. Thus if the argument is close to 1, 
even the roundoff error of the argument causes a substantial relative 
error in the answer, since the function value there is very small. 
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Accuracy : 

I Arguments j Relative Error 

I 1 

I Range | Distribution! RMS | Maximum 

L 

IHELNLO 

I Exclud- 1 I I 

ling I I 1 

|0.5<x I Exponential! 0.0530 | 0,329 

!<2.0 ! ! ! 

L 

IHEL2L0 

I Exclud- ! ! ! 

ling ! ! ! 

|0.5<x ! Exponential! 0.i\^3 j 2,60 
!<2-0 ! ! I 

L 

IHELGLO 

r 

! Exclud- I ! ! 

ling I I ! 

|0,5<x ! Exponential! 0.155 | 0.402 
!<2.0 I ! ! 

r 

I Arguments | Absolute Error 

I ! *10^5 

, ^ 

! Range | Distribution | RMS ! Maximum 

L 

IHELNLO 

r 

!0.5<x ! Uniform ! 0.192 ! 0.507 

|<2.0 I I ! 

IHEL2L0 

!0.5<x I Uniform ! 0.245 | 0.466 
1<2.0 I I ! 

t 

IHELGLO 

r 

|0.5<x ! Uniform ! 0.0318 ! 0.0625 

!<2.0 I ! ! 

L 
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Implementation ; 

• Module size: 360 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the appropriate entry point 
in the table: 



IBM System/ 360 Model Number 



1 Entry 
1 Point 


I 30 1 
1 1 


40 1 


50 


1 65 1 75 


1 IHELNLO 
1 . . ^ 


1 16216 1 

"-+ +- 

1 17315 1 

-+ +- 

1 17284 1 


3926 1 
+- 

4196 1 
+- 

4192 1 


788 


1 178 1 98.5 
.^1 ^ 

1 190 1 107 
1 188 1 105 


1 

IIHEL2L0 
1 ^ ^ 


834 


1 

1 IHELGLO 


828 



Errors Detected 

Error and Exceptional Conditions : 

P: X < (202) 
Local Variables 
None 

Program Interface 
Entry Points 

IHELNLO - Entry for LOG(x) 

P7 = A(PLIST) 
PLIST = A(x) 

A (Target) 

IHELGLO - Entry for LOGlO(x) 

Linkage as for IHELNLO 

IHEL2L0 - Entry for L0G2(x) 

Linkage as for IHELNLO 

Exit Conditions 

Normal. Returns to caller via the link register. 
Abnormal. Branches to EXEP via the entry point IHEERRB. 

Routines Called 

EXEP Execution Error Package 

Global Variables: None 

Comments 

Called by compiled code and IHEATL. 
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TITLE: SHORT FLOAT COMPLEX LOG (IHELNW) 

Program Definitiion 

Purpose and Usage 

Short Float Complex Log is used to calculate the principal value of 
the natural log of z (that is, -pi < imaginary part of result < pi) , 
where z is a short floating-point complex expression. 

Description 

Method: 



Let 



LOG(x ♦ yl) = u + vl 



Then 

U = LOG(ABS(x ♦ yl)) 

= L0G(SQRr(x**2 ♦ y**2)) 
= LOG(x**2 ♦ y**2)/2 

V = ATAN(y,x) 

In computing u, the exponents of x and y are modified if necessary 
to avoid OVERFLOW or UNDERFLOW, with the appropriate correction being 
applied after the logarithm has been taken. 

Effect of an Argument Error : 

Let 



z = r*EXP(hI) 



and 

LOG(z) = u ♦ VI 

Then the absolute error in u is approximately equal to the relative 
error in r* For the absolute error in v (h = ATAN(y,x)), see 
corresponding paragraph for Short Float Real Arctan (IHEATS) . 

Accuracy ; 



Arguments | Relative Error 
I *10« 




I Range | Distribution! RMS 
, ^ + 

I I 

(Expanded | 0.150 

I radially,, | 

I uniform 



I Full 
I Range 
I except 
I within 
|10-« 



I of 



1+0I| 



I 

(round origin | 



Implementation : 

• Module size: 272 bytes 

• Execution times: 
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Approximate execution times in microseconds for the System/360 
models given below are obtained from the table: 



r 

1 1 IBM System/ 360 Model Number 


1 1 30 1 40 1 50 1 65 1 


75 


1 (i) 1 11511 1 3414 1 1078 1 308 | 
, 4. ^ ^ ^ +. 

1 (ii) 1 11688 1 3489 | 1104 | 318 | 
1 (iii) 1 11800 1 3520 | 1117 | 321 | 


183 


190 


193 



(i) ABS(x) and ABS(y) < SQRT(8) *16**31 and either ABS(x) 
or ABS(y) > 16**-30. 

(ii) Either ABS(x) or ABS(y) > SQRT(8)*16**31. 

(iii) ABS(x) and ABS(y) < 16**-30. 
Errors Detected 
Error and Exceptional Conditions ; 

O: X = y = (202) in Short Float Real Log (IHELNS) 
Local Variables 
None 

Program Interface 
Entry Points 

IHELNTO 

P7 = A(PLIST) 
PLIST = A(z) 

A (Target) 

Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

IHELNS Short Float Real Log 
IHEATS Short Float Real Arctan 

Global Variables 

None 

Comments 

Called by compiled code and IHEXXW. 
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TITLE: LONG FLOAT COMPLEX LOG (IHELNZ) 

Program Definition 

Purpose and Usage 

Long Float Complex Log is used to calculate the principal value of 
the natural log of z (that is^ -pi < imaginary part of result < pi) , 
where z is a long floating-point complex expression. 

Description 

Method ; 

Let 



LOG(x + yl) = u + vl 



Then 

u = LOG(ABS(x ♦ yl)) 

= L0G(SQRT(x**2 ♦ y**2)) 
= L0G(x**2 ♦ y**2)/2 

V = ATAN(y,x) 

In computing u, the exponents of x and y are modified if necessary 
to avoid OVERFLOW or UNDERFLOW,^ with the appropriate correction being 
applied after the logarithm has been taken. 

Effect of an Argument Error ; 

Let 

z = r*EXP(hI) 
and 

LOG(z) = u + vl 

Then the absolute error in u is approximately egual to the relative 
error in r. For the absolute error in v (h = ATAN(y,x))t see the 
corresponding paragraph for Long Float Real Arctan (IHEATL) . 

Accuracy ; 
r 



Arguments 



I 



Relative Error 

♦ 101-5 



Range | Distribution j RMS | 



Full I Exponential | 
Range | radially, | 
except I uniform | 
within (round origin) 
I I 



0.0558 



10- « 

of l+OIl 



I 



Maximum 



Implementation ; 

• Module size; 288 bytes 

• Execution times; 
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Approximate execution times in microseconds for System/360 models 
given below are obtained from the table: 



1 IBM System/ 360 Model Number 




1 


30 1 no 1 50 1 65 


1 


75 


(i) 


1 


aUlOl 1 10166 1 2086 | USO 


1 


274 


(ii) 


— +- 

1 
— +- 

1 


^ ^ + 

44398 1 10316 | 2155 | 507 

HHU39 1 10325 | 2156 | 507 


1 
-+- 

1 


290 


(iii) 


290 



(i) 



ABS(x) and ABS(y) < SQRT(8) *16**31 and either ABS(x) 
or ABSCy) > 16**~26. 



(ii) Either ABS(x) or ABS(y) > SQRT(8)*16**31. 

(iii) ABS(x) and ABS(y) < 16**^26. 
Errors Detected 
Error and Exceptional Conditions ; 

O: X = y = (202) in Long Float Real Log (IHELNL) 
Local Variables 
None 

Program Interface 
Entry Points 

IHELNMO 

P7 = A(PLIST) 
PLIST = A(z) 

A (Target) 

Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

IHELNL Long Float Real Log 
IHEATL Long Float Real Arctan 

Global Variables 

None 

Comments 

Called by compiled code and IHEXXZ. 
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TITLE: SHORT FLOAT REAL SIN/COS (IHESNS) 

Program Definition 

Purpose and Usage 

Short Float Real Sin/Cos is used to compute sin(x) and cosCx)^ 
where x is a short floating-point real expression in radians. 

Description 

Method : 

Let 

k = pi/a 
Evaluate 

p = ABS(x)*(l/k) 

where x is in radians using long-precision multiplication to safeguard 
accuracy • Separate p into integer part q and fractional part r, that 
is, p = q ♦ r where < r < 1. 

Define 

ql = q (if SIN is required and x > 0) 

ql = q + 2 (if COS is required) 

ql = q ♦ U (if SIN is required and x < 0) . 

Then for all values of x, each case has been reduced to the computation 
of SIN(k*(ql ♦ r)) = SIN(t), where t > 0. 

Let 

q2 = MOD(ql, 8) 

If q2 = 0, then SIN(t) = SIN(k*r) 

If q2 = 1^ then SIN(t) = COS(k*(l-r)) 

If q2 = 2, then SIN(t) = COS(k*r) 

If q2 = 3, then SIN(t) = SIN(k*(l-r)) 

If q2 = 4, then SIN(t) = -SIN(k*r) 

If q2 = 5, then SIN(t) = -COS(k* (1-r)) 

If q2 = 6, then SIN(t) = -COS(k*r) 

If q2 = 7, then SIN(t) = -SIN(k* (1-r) ) 

Thus it is necessary to compute only SIN(k*rl) or COS(k*rl) where 
rl = r or 1 - r and < rl < !• 

This is performed by using the Chebyshev interpolation polynomials of 
degree 3 in rl**2, with maximum relative error of 2**-28.1 in the sine 
polynomial and 2**-24,6 in the cosine polynomial. 
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Effect of an Argument Error ; 

The absolute error of the answer is approximately equal to the absolute 
error in the argument. Hence, the larger the argument, the larger its 
absolute error and the larger the absolute error of the result. Since 
the function diminishes periodically for both sine and cosine, no 
consistent control of the relative error can be maintained outside the 
range -pi/2 to pi/2 radians. 

Accuracy ; 



i Arguments | 

i i 

, ^. 

I Range | Distribution! 

L 



Absolute Error 
*10« 



RMS 



Maximum 



IHESNSO 



r 

|ABS(x)<| Uniform 


1 0.0557 1 


0.126 




lpi/2 1 








1 + 


— ^ ^«- 





— 


|pi/2< 1 








|ABS(x) 1 Uniform 


1 0.0553 1 


o.ias 




|<10 1 
|10< 1 


--j — j- 










|ABS(x) 1 Uniform 


1 0.0560 1 


o.iaa 




1 <100 1 

L . 









IHECSSO 






r 

1 0<x^pi 
1 

1-1 0<x 

|<0, 

|pi<x 

|<10 

1 

|10< 

|ABS(x) 

|<100 


1 Uniform 
-^ 

1 

1 Uniform 

1 
1 
-+ 

1 

1 Uniform 

1 


1 

~+- 

1 
1 
I 
1 
~+- 
1 
1 
1 


0.0553 1 0.ia9 

1 
0.0571 1 0.15U 

1 

1 
^ 

1 
0.0553 1 0.H»2 

1 




r " 

1 

1 

1 _ ^. 


Arguments 


1 

1 

1 


Relative Error 
*10« 


1 

1 Range 

L 


._ ^ 

1 Distribution | 


RMS 1 Maximum 



IHESNSO 



|ABS(x)<| Uniform 
I pi/2 I 



I 0.198 I 1.40 



I 



I 



Implementation : 

• Module size: 320 bytes 
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• Execution tiroes: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the appropriate entry point 
in the table: 



r 

1 


1 


IBM System/ 360 Model 


Number 




1 

|ABS(x) 
1 


1 30 


1 


HO 


1 


50 


1 


65 


1 


75 


IHESNSO 


|<pi/4 
1 >pi/4 


1 ao9i 

1 t»386 


1 
1 


1120 
1190 


1 

1 


333 

362 


1 

1 


85.0 
92.5 


1 
1 


50.6 
53.8 


IHECSSO 


1 <pi/4 
|>pi/4 

L 


1 ao78 

I 4373 


I 

1 


1115 
118a 


1 
1 


329 
357 


1 

1 


83.6 
91.0 


1 
1 


«»9.9 
53.1 



Errors Detected 

Error and Exceptional Conditions : 

P: IHESNSO, IHECSSO: ABS(x) > 2**1 8*pi (212) 
Local Variables 
None 

Program Interface 
Entry Points 

IHESNSO - Entry for SIN(x) 



P7 


= 


A (PLIST) 


PLIST 


= 


A(x) 

A (Target) 



IHECSSO - Entry for COS(x) 

Linkage as for IHESNSO 

Exit Conditions 

Normal. Returns to caller via the link register. 

Abnormal. Branches to the EXEP via the entry point IHEERRB, 

Routines Called 

EXEP Execution Error Package 

Global Variables: None 

Comments 

Called by compiled code, IHESNW, and IHEEXW. 
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TITLE: LONG FLOAT REAL SIN/COS (IHESNL) 

Program Definition 

Purpose and Usage 

Long Float Real Sin/Cos is used to compute sin(x) and cos(x), 
where x is a long floating-point real expression in radians • 

Description 

Method ; 

Let 

y = ABS(x)/(pi/'*) for x in radians and 

y = q + r^ q integral, < r < 1 
Take 

ql = q for SIN with positive or zero argument 

ql = q + 2 for COS 

ql = q + 4 for SIN with negative argument and 

q2 = MOD(ql, 8) 
Since 

COS(x) = SIN(ABS(x) + pi/2) and 

SIN(-x) = SIN(ABS(x) + pi) 
it is only necessary to find 

SIN (pi/4* (q2 ♦ r)), for < q2 < 7 
Therefore compute: 

SIN(pi/4*r) if q2 = or 4 

COS (pi/a* (1-r)) if q2 = 1 or 5 

COS(pi/4*r) if q2 = 2 or 6 

SIN(pi/4*(l-r)) if q2 = 3 or 7 

SIN(pi/4*rl) where rl is r or (1 - r) is computed by using the Chebyshev 
interpolation polynomial of degree 6 in rl**2, in the range < rl**2 
< 1, with maximum relative error 2**(-58). 

C0S(pi/4*rl) is computed by using the Chebyshev interpolation polynomial 
of degree 7 in rl**2v in the range < rl**2 < 1, with maximum relative 
error 2**(-6U.3). 

Finally, if q2 > 4, a negative sign is given to the result. 

Effect of an Argument Error : 

The absolute error of the answer is approximately equal to the absolute 
error in the argument. Hence, the larger the argument, the larger 
its absolute error and the larger the absolute error of the result. 
Since the function diminishes periodically for both sine and cosine, 
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no consistent control of the relative error can be maintained outside 
the range -pi/2 to pi/2 radians. 

Accuracy : 

IHESNLO 



Arguments 



I Range | Distribution 
I- 



I -pi/2 <x I Uniform 
|<pi/2 I 



Relative Error 



RMS 



Maximum 



0.0542 I 0-381 
I 



IHECSLO 



Arguments 



Range | Distributi on 

-pi/2<x| Uniform 
<pi/2 I 



Relative Error 
♦10^.5 



RMS I Maximum 

0.0604 I 0.168 
I 



Implementation ; 

• Module size: 416 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the appropriate entry point 
in the table: 



IBM System/ 360 Model Number 



I Entry 
I Point 

I 



30 



I 



40 



50 



65 



75 



I IHESNLO I 13654 | 3290 | 661 | 155 | 85.3 
I IHECSLO I 13641 | 3288 | 654 | 153 | 84.3 

L , 

Errors Detected 

Error and Exceptional Conditions : 

P: IHESNLO,^ IHECSLO: ABS(x) > 2**50*pi (203) 
Local Variables 
None 
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Program Interface 

Entry Points 

IHESNLO - Entry for SIN(x) 

P7 = A(PLIST) 
PLIST = A (x) 

A (Target) 

IHECSLO - Entry for COS(x) 

Linkage as for IHESNLO 

Exit Conditions 

Normal. Returns to caller via the link register. 

Abnormal. Branches to the EXEP via the entry point IHEERRC. 

Routines Called 

EXEP Execution Error Package 

Global Variables 

None 

Comments 

Called by compiled code^ IHESNZ, and IHEEXZ. 
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TITLE: SHORT FLOAT REAL HYPERBOLIC SIN/ COS (IHESHS) 

Program Definition 

Purpose and Usage 

Short Float Real Hyperbolic Sin/Cos is used to calculate hyperbolic 
sin(x) or hyperbolic cos(x), where x is a short floating-point real 
expression in radians. 

Description 

Method ; 

For IHESHS, if ABS(x) < 1, use a polynomial approximation of the seventh 
degree. Othencise 

SINH(x) = EXP(x)/2 - 0.5/EXP(x) 

COSH(x) = EXP(x)/2 + 0.5/EXP(x) 

These two versions of EXP(x)/2±0.5/EXP(x) are preferable to the 
equivalent versions of (EXP(x) - l/EXP(x))/2 because^ in floating- 
point, 0.5 has three more significant bits than 1.0. 

Effect of an Argument Error: 

The relative error caused in the result is approximately as follows: 

SINH: The absolute error in the argument divided by TANH(x), that is, 
of the order of the absolute error in the argument for large x, 
or of the relative error in the argument for small x. 

COSH: The absolute error in the argument multiplied by TANH(x), 

that is, of the order of the absolute error in the argument. 

Thus, for large values of x, even the roundoff error of the argument 
causes a substantial relative error in the answer. 

Accuracy : 



I Arguments | 

I I 

I 1. 

I Range | Distribution | 

L 



Relative Error 
*10« 



RMS 



Maximum 



IHESHSO 



|ABS(x)<| Uniform | 0.200 j 0.932 

|1^ I ^ 1^ '. . . 

|1< I II 

|ABS(x)<| Uniform | 0.221 j 0.950 

|2 I II 

L 
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IHECHSO 



JABS(x)<| Onifonn 
|1 I 



I 0.367 
I 



, ^ ^ 

|1< I I 

|ABS(x)<| Uniform | 0.192 

|2 I I 

L 



I 
I 

I 
I 
I 



0.908 



0.700 



Implement-ation ; 

• Module size: 216 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the appropriate entry point 
in the table: 



r 

1 


1 




IBM System/ 360 Model Number 




1 

1 Entry 
1 Point 




30 


1 

1 


UO 


1 

1 


50 


1 65 

1 


1 

1 


75 


1 

1 IHESHSO 

|ABS(X><1 

lABSCx)>l 

1 -. — 




25IH1 
56U7 


1 
1 

1 

1 


692 
1693 


1 

1 

1 

_j 


228 
526 


1 

1 56.3 

1 nm 


1 
1 
1 

^ 1 


33.0 
91.8 


1 

j IHECHSO 

L 




5500 


+ 
1 


1648 


T 
1 


509 


1 139 


t 
1 


88.4 



Errors Detected 

Error and Exceptional Conditions ; 

H: OVERFLOW (300) in Short Float Real EXP (IHEEXS) 
Local Variables: None 
Program Interface 
Entry Points 

IHESHSO - Entry for SINH(x) 

P7 = ACPLIST) 
PLIST = A(x) 

A (Target) 

IHECHSO - Entry for COSH(x) 
Linkage as for IHESHSO 
Exit Conditions: Normal. Retiarns to caller via link register - 
Routines Called 

IHEEXS Short Float Real EXP 
Global Variables: None 
Comments : Called by compiled code. 
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TITLE: LONG FLOAT REAL HYPERBOLIC SIN/COS (IHESHL) 

Program Definition 

Purpose and Usage 

Long Float Real Hyperbolic Sin/Cos is used to calculate hyperbolic 
sin(x) or hyperbolic cos(x), where x is a long floating-point real 
expression in radians. 

Description 

Method ; 

For IHESHL, if ABS(x) < 0.3U65736, compute SINH(x)/x using polynomial 
approximation of degree 5 in x**2, with relative error less than 
2**-61,9. Otherwise t compute 

s = EXPCABSCx)) 
then 

COSH(x) = (s ♦ l/s)/2 

SINH(x) = SIGN(x)*(s - l/s)/2 

Effect of an Argument Error : 

The relative error caused in the result is approximately as follows: 

SINH: The absolute error in the argument divided by TANH(x), that 
is, of the order of the absolute error in the argument for 
large x, or of the relative error in the argument for small 

X. 

COSH: The absolute error in the argument multiplied by TANH(x) , 

that is, of the order of the absolute error in the argument. 

Thus, for large values of x, even the roundoff error of the argument 
causes a substantial relative error in the answer. 

Accuracy ; 



j Argximents | 

I I 

, ^.. 

I Range | Distribution j 
L — . 



Relative Error 
♦10" 



RMS 



Maximum 



IHESHLO 



|ABS(x)<| Uniform 
|0.3U657| 

I 0. 346571 Uniform 

l< I 
|ABS(x)<| 

|5 I 

L 



I 0.0530 I 0.217 
I I 

I 0.0870 I 0.359 



5«* 



IHECHLO 



|ABS(x)<| Uniform | 0.123 | 0,429 | 

|5 I II I 

L J 

Implementation ; 

• Module size: 264 bytes 

• Execution times: 

Approximate execution times in microseconds for System/360 models 
given below are obtained from the appropriate entry point in the 
table: 

J. - ^ 

I I imi System/ 360 Model Number | 
, , 

|ABS(x) I 30 I 40 I 50 I 65 I 75 I 

L J 

IHESHLO 

|<0.347 I 9024 | 2279 | 450 | 101 | 59,0 

|0.347< I I i I I 
|ABS(x) I 18634 | 4338 | 938 | 215 | 125 
|<174.6 i I I I I 

L J 

IHECHLO 

r 1 

|<174,6 I 18493 | 4300 | 924 | 211 | 123 | 

L J 

Errors Detected 

Error and Exceptional Conditions ; 

H: OVERFLOW (300) in Long Float Real EXP (IHEEXL) 
Local Variables 
None 

Program Interface 
Entry Points 

IHESHLO - Entry for SINH(x) 

P7 = A(PLIST) 
PLIST = A(x) 

A (Target) 

IHECHLO - Entry for COSH(x) 

Linkage as for IHESHLO 
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Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

IHEEXL Long Float Real EXP 
Global Variables 
None 
Comments 
Called by compiled code. 
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TITLE: SHORT FLOAT COMPLEX SIN/COS (IHESNW) 

Program Definiiiion 

Purpose and Usage 

Short Float Complex Sin/Cos is used to calculate sin(z)t hyperbolic 
sin(z), cos(z), and hyperbolic cos(z),, where z is a short floating- 
point complex expression in radians. 

Description 

Method: 

Let 

z = X + yl 
Then 

REAL (SIN (z)) = SIN(x)*COSH(y) and 

IMAG(SIN(z)) = COS(x)*SINH(y) 

REAL(COS(z)) = COS(x)*COSH(y) and 

IMAG(COS(z)) = -SIN(x)*SINH(y) 

REAL(SINH(z)) = COS(y) *SINH(x) and 

IMAG(SINH(z)) = SIN(y)*COSH(x) 

REAL (COSH (z)) = COS(y) *COSH (x) and 

IMAG(COSH(z)) = SIN(y)*SINH(x) 

To avoid making calls to evaluate SINK and COSH separately, and thus 
frequently having to evaluate EXP twice for the same argument, SINH(u) 
is computed as follows: 

1. u > 0,3465736 

SINH(u) = (EXP(U) - l/EXP(u))/2 

2. < u < 0.3465736 

SINH(u)/u is approximated by a polynomial of the form 
a ♦ al*u**2 + a2*u**4 (which has a relative error of 
less than 2**-26.4). 

3. u < 

SINH(u) = ~SINH(-u) 

Then 

COSH(u) = SINH(ABS(u)) ♦ 1/EXP(ABS (u) ) 

Effect of an Argument Error : 

Combine the effects on SIN, COS, SINH, and COSH according to the method 
of evaluation described in the above paragraph. 
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Accuracy ; 

I Argiiments | Relative Error | 

I I *10*^ I 

I ^ , 

I Range | Distribution | RMS | Maximum | 

L J 

IHESNTO 

J. ^ 

|ABS(x) I Uniform | 0.721 | 2.07 j 

1^10. Ill I 

|ABS(y)<| I I I 

|1 I i I I 

L J 

IHESHTO 

|ABS(x) I Uniform | 0.561 | 1.86 | 

1^10, III I 

|ABS(y)<| I I I 

|1 I II I 

L J 

IHECSTO 

J. ^ 

|ABS(x) I Uniform | 0.5^6 | 2.00 | 

1^10. Ill I 

|ABS(y) III I 

|<20 III I 

L J 

IHECHTO 

|ABS(x) I Uniform | 0.558 | 2.35 | 

1^10. Ill I 

|ABS(y) III I 

|<20 III I 

L J 

Implementation : 

• Module size: 320 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the appropriate entry point 
in the tables: 
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SIN^COS : ABS(y)>0. 3465736 
SINH^COSH : ABS(x)>0. 3465736 



r 

1 1 IBM System/ 360 Model 


Number 


^ 


1 Entry | 30 | «»0 | 50 | 

1 Point 1 1 1 1 
, ^ ^ ^ ^_. 

IIHESNTO 1 15826 | 1508 | 648 | 

IIHECSTO 1 15898 | 1518 | 653 | 

IIHESHTO 1 15930 j 1520 | 653 | 
, ^ 1 ^ ^.. 

IIHECHTO 1 15900 | 1519 | 655 | 


65 1 

I 
+- 

363 1 
+- 

366 1 
+- 

366 1 
+- 

367 1 


75 1 


75 1 
+ 

225 1 

226 1 


227 1 



SIN, COS : ABS(y)<0. 3465736 
SINH^COSH : ABS(x)<0. 3465736 



j IHESNTO 

I 

I IHECSTO 

I IHESHTO 

I IHECHTO 

L 



I 16896 I 1585 | 674 | 381 | 232 

I 16968 I 1595 | 679 | 384 | 234 

I 17000 I 1596 I 679 | 384 | 234 

I 16970 I 1595 I 681 | 384 | 235 



Errors Detected 

Error and Exceptional Conditions ; 

0: IHESNWS, IHECSTO: ABS(x) > 2**18*pi (212) in Short Float Real 

Sin/Cos (IHESNS) 
H: OVERFLOW (300) in Short Float Real EXP (IHEEXS) 

Local Variables: None 

Program Interface 

Entry Points 

IHESNTO - Entry for SIN(z) 

P7 = A(PLIST) 
PLIST = A (Target) 

IHECSTO - Entry for COS(z) 

Linkage as for IHESNTO 

Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 



IHESNS Short Float Real Sin/Cos 
IHEEXS Short Float Real EXP 

Global Variables: None 

Comments : Called by compiled code. 



59 



TITLE: LONG FLOAT COMPLEX SIN/COS (IHESNZ) 

Program Definition 

Purpose and Usage 

Long Float Complex Sin/Cos is used to calculate sin(z), hyperbolic 
sin(z)^ cos(z)^ or hyperbolic cos(z), where z is a long floating-point 
complex expression in radians. 

Description 

Method: 

Let 

z = X + yl 
Then 

REAL ( SIN (z)) = SIN (x) ♦COSH (y) and 

IMAG(SIN(z)) = COS (x) ♦SINK (y) 

REAL(COS(z)) = COS(x)*COSH(y) and 

IMAG(COS(z)) = -SIN(x)*SINH(y) 

REAL (SINE (z)) = COS (y) *SINH(x) and 

IMAG(SINH(z)) = SIN (y) ♦COSH (x) 

REAL (COSH (z)) = COS (y) ♦COSH (x) and 

IMAG(COSH(z)) = SIN(y)^SINH(x) 

To avoid making calls to evaluate SINH and COSH separately,, and thus 
frequently having to evaluate EXP twice for the same argument, SINH(u) 
is computed as follows: 

1. u > 0.3465736 

SINH(u) = (EXP(u) - l/EXP(u))/2 

2. < u < 0.3465736 

SINH(u)/u is approximated by a polynomial of the fifth degree 
in u^^2 which has a relative error of less than 2^^-61.8. 

3. u < 

SINH(u) = -SINH(~u) 

Then 

COSH(u) = SINHCABSCu)) + 1/EXP(ABS (u) ) 

Effect of an Argument Error : 

Combine the effects on SIN, COS, SINH, and COSH according to the method 
of evaluation described in the above paragraph. 
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Accuracy : 

J. ^ 

I Arguments | Relative Error | 

, ^ 1 

I Range | Distribution | RMS | Maximum | 

L . . . J 

IHESNMO 

^ ^ 

|ABS(x) I Uniform | 2.11 J 82.4 | 
1^10. Ill I 

|ABS(y)<| I I I 

|1 I I i i 

IHESHMO 

|ABS(x) I Uniform | 0.180 | 2.31 | 

1^10, I II I 

|ABS(y)<I I I I 

|1 I II I 

L J 

IHECSMO 

JABSCx) I Uniform | 0.389 | 6.33 | 

1^10. I 11 I 

|ABS(y)<| I I I 

|1 I II I 

L . J 

IHECHMO 

r 1 

|ABS(x) I Uniform | 1.11 | 19.4 | 

1^10, I II I 

|ABS(y) I II I 

|<20 I II I 

L J 

Implementation ; 

• Module size: 368 bytes 

• Execution tiroes: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the appropriate entry point in 
the tables: 
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SIN^COS : ABS(y)>0, 3465736 
SINH,COSH : ABS(x)>0. 3465736 



r 
1 

1 -...•«.. ^ 


1 


IBM Systein/360 Model 


Number 


1 


1 

1 Entry 
1 Point 


1 
1 

1 

1 
1 
1 


30 1 HO 1 50 1 
1 1 1 

«»658i» 1 11267 1 2363 | 

46656 1 11294 | 2373 I 
+ ^ +-. 

46726 1 11317 | 2378 | 
^ ^ 4... 

U6640 1 11304 1 2374 | 


65 1 

1 
+- 

552 1 
+- 

555 1 
+- 

557 1 
+- 

556 1 


75 1 


1 IHESNMO 

1 -«—.«.--. -.— . 


313 1 
+ 

315 1 


1 

1 IHECSMO 

1^ — «-. « 


1 

1 IHESHMO 

1 


316 1 


I 

1 IHECHMO 

L 


316 1 

J 



SIN, COS : ABS(y)<0. 3465736 
SINH^COSH : ABS(x)<0. 3465736 

r T 

I IHESNMO I 54173 | 13141 | 2656 | 612 | 345 

(IHECSMO I 54245 | 13168 | 2666 | 615 | 347 

I IHESHMO I 54325 | 13191 | 2671 | 617 | 347 

I IHECHMO I 54247 | 13177 | 2667 | 616 | 348 

L J 

Errors Detected 

Error and Exceptional Conditions ; 

0: IHESNZS, IHECSMO: ABS(x) > 2**50*pi (203) in Long Float Real 

Sin/Cos (IHESNL) 
H: OVERFLOW (300) in Long Float Real EXP (IHEEXL) 

Local Variables: None 

Program Interface 

Entry Points 

IHESNMO - Entry for SIN(z) 

P7 = A(PLIST) 
PLIST = A(z) 

A (Target) 

IHECSMO - Entry for COS(z) 

Linkage as for IHESNMO 

Exit Conditions: Nonnal. Returns to caller via the link register. 

Routines Called 

IHESNL Long Float Real Sin/Cos 
IHEEXL Long Float Real EXP 

Global Variables: None 

Comments : Called by compiled code. 
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TITLE: SHOPT FLOAT REAL SQRT (IHESQS) 

Program Definition 

Purpose and Usage 

Short Float Real SQRT is used to compute the square root of x, where 
X is a short floating-point real expression. 

Description 

Method: 

If 

X = 0,, SQRT(x) = 
Otherwise^ let 

X = 16**(2*p ♦ q)*f 
where p is an integer^ q = or 1, and 1/16 < f < 1. Then 

SQRT(x) = 16**(p ♦ q)*z 
where 

z = SQRT(f) if q = 

z = SQRT(f)/4 if q = 1 

An initial approximation, yO, is taken in the hyperbolic form 

a ♦ b/(c ••• f) with different sets of constants for the two cases: 

1. q = a = 1,80713 

b = 1.57727 

c = 0.95m82 

The maximum relative error in this range is less than 
2** (-5. 4U) g with an exact fit at f = 1 to guard as far as 
possible against loss of the last hexadecimal digit v^en 
f is nearly 1, 

2. q = 1 a = 0,428795 

b = -0.3430368 

c = 0.877552 

The maximum relative error in this range is less than 
2**(-6)*f**(-l/8). 

Then 

yl - 16** (p ♦ q)*yO 
Two Newton-Raphson iterations then yield: 

y2 = (yl + x/yl)/2 

SQRT(x) = y2 * (x/y2 - y2)/2 
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For case q = 0, the final relative error from this algorithm is less 
than 2**(-24.7), and, for case q = 1, less than 2**(-29), 

Effect of an Argument Error : 

The relative error caused in the result is approximately half the 
relative error in the argument. 

Accuracy ; 



Arguments 



I 
I 



Relative Error 
*10« 



RMS 



I Range | Distribution j 

I ^ ^ ^., 

I Full I Exponential! 0.230 | 
I Range | | | 

L , ., 



Maximum 
0.924 



Implementation ; 

• Module size: 168 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the table: 



1 




IBM Systein/360 Model Number 




1 


30 


1 40 1 50 1 65 1 


75 


1 


31«»0 


1 793 1 227 1 68.4 | 


40.7 



Errors Detected 

Error and Exceptional Conditions ; 

P: X < (200) 
Local Variables 
None 

Program Interface 
Entry Points 

IHESQSO 

P7 = A(PLIST) 
PLIST = A(x) 

A (Target) 

Exit Conditions 

Normal. Returns to caller via the link register. 
Abnormal. Branches to the EXEP via IHEERRB. 

Routines Called 

EXEP Execution Error Package 
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Global Variables 

None 

Comments 

Called by compiled code^ IHEABW^ and IHESQW. 
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TITLE: LONG FLOAT REAL SQRT (IHESQL) 

Program Definition 

Purpose and Usage 

Long Float: Real SQRT is used to compute the square root of x, where 
X is a long floating-point real expression. 

Description 

Method ; 

If 

X = 0^ SQRT(x) = 
otherwise t let 

X = 16**(2*p - q)*f 
where p is an integer, 

q = or 1, and 1/16 < f < 1 
Then 

SQRT(x) = 16**p*2**(-2*q)*SQRT(f) 

An initial approximation,, yO, is taken by using (2/9 ♦ 8/9*f) for 
SQRT(f), Multiplication by 2** (-2) when q = 1 is accomplished by using 
the HALVE instruction twice. The maximum relative error of this 
approximation is 1/9. 

Four Newton- Raphson iterations of the form y(n+l) = (yn + x/yn)/2 are 
then applied, two in short precision and two in long precision, the 
last being computed as 

SQRT(x) = y3 ♦ (x/y3 - y3)/2 

to minimize the truncation error. 

The maximum relative error in the result from this algorithm is 
2**(-65.7). 

Effect of an Argument Error ; 

The relative error caused in the result is approximately half of the 
relative error in the argument. 

Accuracy ; 

^ ^ 

Argtiments ] Relative Error 
I *10^s 
^ 

Range | Distribution! RMS | Maximum 
|x>10-s2| Exponential! 0.0276 | 0.12U 
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Implementation : 

• Module size: 160 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the table: 



1 




IBM Syst,ein/360 Model Number 




1 


30 


1 40 1 50 1 65 1 


75 


1 


8282 


1 1733 1 376 1 97.7 | 


57.2 



Errors Detected 

Error and Exceptional Conditions : 

P: X < (200) 
Local Variables 
None 

Program Interface 
Entry Points 

IHESQLO 

P7 = A(PLIST) 
PLIST = A(x) 

A (Target) 

Exit Conditions 

Normal. Returns to caller via the link register. 
Abnormal. Branches to EXEP via IHEERRB. 

Routines Called 

EXEP Execution Error Package 

Global Variables 

None 

Comments 

Called by compiled code, IHEABZ, and IHESQZ. 
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TITLE: SHORT FLOAT COMPLEX SQRT (IHESQW) 

Program Definitdon 

Purpose and Usage 

Short Float Complex SQRT is used to calculate the principal value of 
the square root of 2 (that is^ -pi/2 < argument of result < pi/2) , 
where z is a short f loating-^point complex express ion • 

Description 

Method ; 

Let 

z = X + yl and 

SQRT(z) = u ♦ VI 
If 

1. X = y = 
Then 

u = V = 

2. X > 
Then 

u = SQRT((ABS(x) + ABS(x + yI))/2) and 

V = y/(2*u) 

3. X < 
Then 

u = y/(2*v) and 

V = S(y)*SQRT((ABS(x) + ABS (x + yI))/2) 

where 

S(y) = 1 if y > 
= -1 if y < 

Effect of an Argument Error ; 

Let 

2 = r*EXP(hI) and 

SQRT(z) = s*EXP(kI) 

Then the relative error in s is approximately half the relative error 
in r, and the relative error in k is approximately equal to the relative 
error in h. 
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Accuracy ; 



Arguments | Relative Error 
I *10« 



Range (Distribution! RMS | 



Full I Exponential! 0,513 

Range | radially^ | 

! uniform | 

! round | 

! origin \ 



Maximum 
1,51 



Implementation : 

• Module size: 152 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the table: 



1 


IBM Systein/360 Model Number 




1 


30 1 40 1 50 1 65 1 


75 


1 


11130 1 3023 1 1006 | 265 ] 


16«t 



Errors Detected 

Error and Exceptional Conditions : 

I: OVERFLOW (300) 

H: OVERFLOW (300) in Short Float Complex ABS (IHEABW) 

Local Variables: None 

Program Interface 

Entry Points 

IHESQTO 

P7 = A(PLIST) 
PLIST = A(z) 

A (Target) 

Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

IHEABW Short Float Complex ABS 
IHESQS Short Float Real SQRT 

Global Variables: None 

Comments 

Called by compiled code. 
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TITLE: LONG FLOAT COMPLEX SQRT (IHESQZ) 

Proqraro Definition 

Purpose and Usage 

Long Float Complex SQRT is used to calculate the principal value of 
the square root of z (that is, -pi/2 < argument of result < pi/ 2) , 
where z is a long floating-point complex expression. 

Description 

Method ; 

Let 

z = X + yl and 

SQRT(z) = u ♦ VI 

If 

!• X = y = 

Then 

u = V = 

2. X > 
Then 

u = SQRT((ABS(x) + ABS(x -i- yI))/2) and 

V = y/(2*u) 

3. X < 
Then 

u = y/(2*v) and 

V = S (y) ♦SQRT ( CABS (x) + ABS (x + yI))/2) 

where 

S(y) = 1 if y > 
= -1 if y < 

Effect of an Argument Error ; 

Let 

z = r*EXP(hI) and 

SQRT(z) == s*EXP(kI) 

Then the relative error in s is approximately half the relative error 
in r,^ and the relative error in k is approximately equal to the relative 
error in h. 
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Accuracy ; 



Arguments | Relative Error 
I *10" 



Range | Distribution! RMS | Maximum 



Full 
Range 



Exponential j 0. 263 
radially, | 



uniform 

round 

origin 



l,5t* 



Implementation ; 

• Module size: 14 ^ bytes 

• Execution times: 

Approximate execution times in microseconds for System/360 models 
given below are obtained from the table: 



r 


1 


IBM Systein/360 Model Number 






1 


30 1 HO 1 50 1 65 1 


75 


L 


1 


26996 1 5957 | 1352 | 3tH | 


203 



Errors Detected 

Error and Exceptional Conditions : 

I: OVERFLOW (300) 

H: OVERFLOW (300) in Long Float Complex ABS (IHEABZ) 

Local Variables: None 

Program Interface 

Entry Points 

IHESQMO 

P7 == A(PLIST) 
PLIST = A(z) 

A (Target) 

Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

IHEABZ Long Float Complex ABS 
IHESQL Long Float Real SQRT 

Global variables: None 

Comments 

Called by compiled code. 
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TITLE: SHORT FLOAT REAL TAN (IHETNS) 

Program Definition 

Purpose and Usage 

Short Float Real Tan is used to compute tan(x), where x is a short 
floating-point real expression in radians. 

Description 

Method ; 

Evaluate 

p = (4/pi)*ABS(x) 

Let q and r be respectively the integral and fractional parts of p. 

If q is even, then s = r 

If q is odd, then s = 1 - r 

Let 

ql = M0D(q,4) 
Then 

If ql = 0, then TAN(ABS(x)) = TAN(pi*s/4) 

If ql = 1, then TAN CABS (x)) = COT(pi*s/4) 

If ql = 2, then TAN(ABS(x)) = -COT(pi*s/4) 

If ql = 3, then TANCABSCx)) = -TAN(pi*s/«») 
Compute TAN(pi*s/4) and COT(pi*s/4) as the ratio of two polynomials: 

TAN(pi*s/U) = s*p(s**2)/q(s**2) 

COT(pi*S/a) = q(s**2)/(s*p(s**2)) 
where 

p(s**2) = 212.58037 - 12. 559912*s**2 

q(s**2) = 270.665736 - 71. 6U5273*s**2 + s**4 
Finally, if x < 0, 

TANCx) = -TANCABSCx)) 

Effect of an Argument Error ; 

The absolute error of the answer is approximately equal to the absolute 

error of the argument multiplied by CI + TANCx) **2). Hence if x is 

near an odd multiple of pi/2, an argument error will produce a large 
absolute error in the answer. 

The relative error in the result is approximately equal to twice the 
absolute error in the argument divided by SINC2*x). Hence, if x is 
near a multiple of pi/2, an argument error will produce a large relative 
error in the result. 
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Accuracy ; 



r 
1 
1 


Arguments 


1 

1 


1 

Relative Error | 

*10e 1 


1 T 

1 Range | Distribution] 
|ABS(x)<| Uniform | 


RMS 
0.319 


1 

-+- 


1 

Maximum | 
1,92 1 


lpi/4 


1 

-+- 




1 




-+- 




1 pi/a< 
|ABS(x) 


Uniform 


0.U65 


1.24 1 


|<1.5 














1 

|pi/4< 
|ABS(x) 


-+- 


Uniform 


--+- 


3.ia 


-+- 


170* I 


|<pi/2 

1 -^— «.. 


1 




._j 




— X^. 




|pi/2< 
|ABS(x) 


t 


Uniform 


t 


1.25 


T 


70.6* 1 


|<10 














1 

|10< 


-+- 


.——«————- 


--+- 




-+- 




1 ABS (x) 




Uniform 




3.57 




205* 1 


|<100 















♦These maximum errors are those encoiantered in a sample of 5000 
points; each figure depends very much on the particular points 
encountered near the singularities of the function. 

Implementation ; 

• Module size: 280 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the table: 



I 



IBM System/360 Model Number 



\ 

|ABS(x) I 30 I UO I 50 I 65 I 75 

|<pi/4 I 4429 I 1172 I 336 | 85.8 | 51.0 
|>pi/4 I 4788 I 1262 | 368 | 95.1 | 55.0 

L 

Errors Detected 

Error and Exceptional Conditions ; 

P: ABS(x) > 2**18*pi (213) 
I; OVERFLOW (300) 

Local Variables 

None 
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Procrram Interface 
Entry Points 
IHETNSO 



P7 


= 


A (PLIST) 


PLIST 


~ 


A(x) 

A (Target) 



Exit Conditions 

Normal • Returns to caller via the link register. 
Abnormal. Branches to EXEP via the entry point IHEERRB. 

Routines Called 

EXEP Execution Error Package 

Global Variables 

None 

Comments 

Called by compiled code and IHETNW. 
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TITLE: LONG FLOAT REAL TAN (IHETNL) 

Program Definition 

Purpose and Usage 

Long Float Real Tan is used to compute tan(x), where x is a long 
floating-point real expression in radians. 

Description 

Method ; 

Evaluate 

p = (4/pi)*ABS(x) 
Let q and r be respectively the integral and fractional parts of p. 

If q is even^ then s = r. 

If q is oddr then s = 1 - r. 
Let 

q2 = MOD(q,U) 
Then 

If q2 = 0, then TAN(ABS(x)) = TAN(pi*s/4) 

If q2 = 1, then TAN(ABS(x)) = COT(pi*s/a) 

If q2 = 2, then TAN(ABS(x)) - -C0T(pi*s/4) 

If q2 = 3, then TAN(ABS(x)) = -TAN (pi* s/U ) 
Compute TAN(pi*s/U) and COT(pi*s/4) as the ratio of two polynomials: 

TAN(pi*s/U) = s*p(s**2)/q(s**2) 

COT(pi*s/U) = q(s**2)/(s*p(s**2)) 

where p(s**2) is of degree 3 and q(s**2) is of degree U in s**2t and 
maximum relative error is 3.«**10**-19, 

Finally^ if x < 0,, 

TAN(x) = -TAN(ABS(x)) 

Effect of an Argument Error : 

The absolute error in the result is approximately equal to the absolute 
error in the argument multiplied by (1 + TAN(x)**2). Hence, if x is 
near an odd multiple of pi/2, an argument error will produce a large 
absolute error in the result- 

The relative error in the result is approximately equal to twice the 
absolute error in the argument divided by SIN(2*x). Hence, if x is 
near a multiple of pi/2, an argument error will produce a large relative 
error in the result. 
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Accuracy ; 
r 



1 Arguments 
1 

1 — — _ _ — «- — . 


1 Relative Error 
1 *10^s 


1 T 

1 Range | Distribution! RMS | 


Maximum 


1 T 

|ABS(x)<| Uniform 


T T 
1 0.091 1 


0.530 


Ipi/i* 1 

1 ^ 

Ipi/iK 1 

|ABS(x) 1 Uniform 


1 0.i\31 1 


2.31 


|<1.5 1 

, ,, 

|pi/4< 1 

|ABS(x) 1 Uniform 


1 7.75 1 




U16* 


|<pi/2 1 

1 1 

|pi/2< 1 
ABS(x) 1 Uniform 


1 18.3 1 


1140* 


— 1 — 

o 1 V 

VI 1 H 


--j !-- 






|ABS(x) 1 Uniform 


1 271 1 


13400* 


|<100 1 







*These maximum errors are those encoimtered in a sample of 5000 
points; each figure depends very much on the particular points 
encountered near the singularities of the function. 

Implementation : 

• Module size: 352 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained frcan the table: 



|ABS(x) 
I- 



1 



30 



IBM System/360 Model Number 
I 40 I 50 I 65 I 



75 



|<pi/4 I 15440 I 3622 | 687 | 154 | 87.2 
|>pi/4 I 16130 I 3817 | 747 | 169 | 93.9 

L 

Errors Detected 

Error and Exceptional Conditions : 

P: ABS(x) > 2**50*pi (204) 
I: OVERFLOW (300) 

Local Variables 

None 
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Program Interface 

Entry Points 

IHETNLO 

P7 = A(PLIST) 
PLIST = A(x) 

A (Target) 

Exit Conditions 

Normal, Returns to caller via the link register. 
Abnormal. Branches to EXEP via the entry point IHEERRB. 

Routines Called 

EXEP Execution Error Package 

Global Variables 

None 

Comments 

Called by compiled code and IHETNZ. 
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TITLE: SHORT FLOAT REAL HYPERBOLIC TAN (IHETHS) 

Program Description 

Purpose and Usage 

Short Float Real Hyperbolic Tan is used to calculate hyperbolic tan(x), 
where x is a short floating-point real expression in radians. 

Description 

Method: 

1. ABS(x) < 2**-12 
Return x as result. 

2. 2**-12 < ABS(x) < 0.54931 

Use a transformed continued fraction of the form: 

TANH(x)/x = 1 - ((x**2 + a)/(x**2 + b + c/x**2)) 
with relative error less than 2**-27. 

3. 0.54931 < X < 9.011 

TANH(x) = 1 - 2/(EXP(2*x) + 1) 

4. X > 9.011 
Return result 1. 

5. X < -0.54931 

TANH(x) = -TANH(-X) 

Effect of an Argument Error: 

The relative error caused in the result is approximately twice the 
absolute error in the argument divided by SINH(2*x). Thus for small 
values of x, it is of the order of the relative error in the argument, 
and as x increases, the effect of the argument error is diminished. 

Accuracy : 

r 1 
Arguments | Relative Error 
I *10e 
^ 

Range | Distribution! RMS | Maximum 

|-0.5<x I Uniform | 0.174 | 0.867 
l<0^5 J^^ I I 

|-9<x<9 I Uniform | 0.0720 | 0.782 

L J 

Implementation : 

• Module size: 200 bytes 

• Execution times: 
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Approximate execution times in microseconds for the System/360 
models given below are obtained from the table: 



r 

1 1 


IBM System/360 Model Number 


1 1 30 

1 


1 40 1 50 1 65 1 75 


1 

|ABS(x)< 1 791 

|2-" 1 

^ ^ 

\2-^^< 1 3033 
|ABS(x)<0.5| 


1 263 1 102 1 28.7 | 21-7 
1 785 1 231 1 64.1 1 U3.9 


|0.5<x 1 5934 
|<9 1 


1 1805 1 562 1 152 | 117 



|3^9 I 1095 I 363 1 139 t '♦0.5 | 35.2 

L , 



Errors Detected 

None 

Local Variables 

None 

Program Interface 

Entry Points 

IHETHSO 

P7 
PLIST 



A(PLIST) 
A(x) 

A (Target) 



Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

IHEEXS Short Float Real EXP 
Global Variables 
None 
Comments 
Called by compiled code and IHETNW. 
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TITLE: LONG FLOAT REAL HYPERBOLIC TAN (IHETHL) 

Program Definition 

Purpose and Usage 

Long Float Real Hyperbolic Tan is used to calculate hyperbolic tan(x), 
where x is a long floating-point real expression in radians. 

Description 

Method: 

ABS(x) < 0,54931 



1. 



2. 



Compute TANH(x)/x using a rational approximation^ with relative 
error less than 2**-6U.5« 



0.51*931 < X < 20.101 

TANH(x) = 1 - 2/(EXP{2*x) ♦ 1) 

3. X > 20.101 
Return result 1. 

4. X < -0.54931 

TANH(x) = -TANH(-x) 

Effect of an Argument Error : 

The relative error caused in the result is approximately twice the 
absolute error in the argument divided by SINH(2*x). Thus for small 
values of X, it is of the order of the relative error in the argument, 
and as x increases, the effect of the argument error is diminished. 

Accuracy ; 



Arguments 



I 
I 



Relative Error 
♦10^5 



Range | Distribution j RMS | Maximum 



ABS(x)<| Uniform 

0.549311 

1 

0.549311 
<ABS (x) I 
<5 I Uniform 



0.0440 I 0.211 



I 
I 

I I 

I I 

I 0.0250 I 



0.199 



Implementation ; 

• Module size: 280 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the table: 



80 



r ~ * ■ 

I I IBM System/360 Model Niimber 

, 

|ABS(x) I 30 I aO I 50 I 65 I 75 

|<0,5t*9 I 12745 I 3030 | 56U | 123 | 67.9 

|0-549< I i I I I 
|ABS(x) I 16U00 I 3918 | 878 | 205 | 119 
l<20^1^^ I I 1^ I I 

|>20.1 I 1239 I 372 | 135 | 39,3 | 25.5 

L , 

Errors Detected 

None 

Local Variables 

None 

Program Interface 

Entry Points 

IHETHLO 

P7 = ACPLIST) 
PLIST = A(x) 

A (Target) 

Exit conditions 

Normal. Returns to caller via the link register. 

Routines Called 

IHEEXL Long Float Real EXP 

Global Variables 

None 

Comments 

Called by compiled code and IHETNZ. 
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TITLE: SHORT FLOAT COMPLEX TAN/HYPERBOLIC TAN (IHETNW) 

Program Definition 

Purpose and Usage 

Short Float Complex Tan/Hyperbolic Tan is used to calculate tan(z) 
or hyperbolic tan(z)» where z is a short floating-point complex 
expression in radians. 

Description 

M ethod; 

Let 

z = X + yl 
Then 

REAL(TAN(z)) = TAN(x)*(l - TANH(y) **2)/(l + (TAN(x)*TANH(y) )**2) 

IMAG(TAN(Z)) = TANH(y)*(l ♦ TAN(x)**2)/(l ♦ (TAN(x)*TANH(y) )**2) 

TANH(z) = -(TAN(zI))I 

Effect of an Argument Error: 

The absolute error caused in the result is approximately equal to the 
absolute error in the argument divided by ABS(COS(z)**2) for IHETNTO 
or divided by ABS(COSH(z)**2) for IHETHTO. The relative error caused 
in the result is approximately twice the absolute error in the argument 
divided by ABS (SIN(2*z)) for IHETNTO, or divided by ABS(SINH(2*z) ) 
for IHETHTO. 

Accuracy : 

I Arguments | Relative Error | 

I I *io« I 

, ^ , 

I Range | Distribution | RMS | Maximum | 
L J 

IHETNTO 

r 1 

|ABS(x)<| Uniform | 0.^30 | 1.67 | 

|1 I i I I 

|ABS(y)<| I I I 

|9 I II I 

L J 

IHETHTO 

r T 

|ABS(x)<| Uniform | 0.430 I 1.45 | 

|9 I II I 

lABSCyXI I I I 

|1 I II I 

L J 

Implementation : 

• Module size: 184 bytes 
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• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the formulas: 

IHETNTO: a ••• time for execution of IHETHS with argument y 
IHETHTO: b ♦ time for execution of IHETHS with argument x 



r 


1 




IBM Systein/360 Model Number 






1 


30 


1 no 1 50 1 65 1 


75 


1 a 

t b 

L 


1 

1 


9094 
9197 


1 2310 1 696 1 186 | 
1 245«t 1 716 1 191 1 


111 
115 



Errors Detected 

Error and Exceptional Conditions : 

I: OVERFLOW (300) 

O: ABS(u) > 2**18*pi (213) 

where u = x for IHETNTO 
u = y for IHETHTO 

H: OVERFLOW (300) in Short Float Real Tan (IHETNS) 
Local Variables 
None 

Program Interface 
Entry Points 

IHETNTO - Entry for TAN(z) 



P7 


= 


A (PLIST) 


PLIST 


= 


A(z) 

A (Target) 



IHETHTO - Entry for TANH(z) 

Linkage as for IHETNTO 

Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

IHETHS Short Float Real Hyperbolic Tan 
IHETNS Short Float Real Tan 

Global Variables 

None 

Comments 

Called by compiled code. 
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TITLE: LONG FLOAT COMPLEX TAN/HYPERBOLIC TMi (IHETNZ) 

Program Definition 

Purpose and Usage 

Long Float Complex Tan/Hyperbolic Tan is used to calculate tan(z) or 
hyperbolic tanCz), where z is a long floating-point complex expression 
in radians • 

Description 

Method ; 

Let 

z = X + yl 
Then 

i^EAL(TAN(z)) = TAN(X)*(1 - TANH(y) **2)/(l ♦ (TAN(x)*TANH(y) ) **2) 

IMAG(TAN(2)) = TANH(y)*(l ♦ TAN(x)**2)/(l + (TAN(x)*TANH(y) )**2) 

TANH(z) = -CTANCzI))! 
Accuracy ; 

I Arguments | Relative Error | 
I I *10^s I 

I ^ , 

I Range | Distribution | RMS | Maximum | 

L ^ J 

IHETNMO 

|ABS(x)<| Uniform | 0.139 | 1,11 | 

|1 I II I 

|ABS(y)<| j I I 

1^ I II I 

L J 

IHETHMO 

J. * ^ 

|ABS(x)<| Uniform 1 0.137 | 0.980 | 

|9 I 11 I 

|ABS(y)<| I I I 

|1 I II I 

L ^* J 

Implementation ; 

• Module size; IS** bytes 

• Execution times; 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the appropriate entry point 
in the tables; 
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r 

1 


1 




IBM System/ 360 Model Number 




1 

1 ABS(y) 


1 


30 


1 «»0 1 50 1 65 1 


75 


IHETNMO 



^ ^ 

|<0.549 I 40843 I 9625 | 1844 ) 411 | 233 | 

|0.549< I I I i I I 
|ABS(y) I 44498 | 10513 | 10365 | 493 | 284 | 
l<20.1 1^ II II I 

|>20.1 I 29337 I 6967 | 1415 | 327 | 190 | 

.J 

I 
I 

I 

J 





r ~ 
1 

1 — -.— 


1 


IBM System/360 Model Number 




1 

1 ABS(x) 


1 


30 1 40 1 50 1 65 1 


75 


IHETHMO 


|<0.549 
1 


1 


41122 1 9709 1 1871 | ftl9 | 


236 



|0.549< I I I I I 
|ABS(x) I 44777 | 10597 | 2185 | 501 | 287 
l<20^1 I I I I I 

|>20,1 I 29616 I 7051 | 1442 | 334 | 193 

L 

Errors Detected 

Error and Exceptional Conditions ; 

I: OVERFLOW (300) 

O: ABS(u) > 2**50*pi (204) 

where u = x for IHETNMO 
u = y for IHETHMO 

H: OVERFLOW (300) in Long Float Real Tan (IHETNL) 
Local Variables 
None 

Program Interface 
Entry Points 

IHETNMO - Entry for TAN(z) 



P7 


= 


A (PLIST) 


PLIST 


= 


A(z) 

A (Target) 



IHETHMO - Entry for TANH(z) 
Linkage as for IHETNMO 
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Exit Conditions 

Normal, Returns to caller via the link register. 

Routines Called 

IHETHL Long Float Real Hyperbolic Tan 
IHETNL Long Float Real Tan 

Global Variables 

None 

Comments 

Called by compiled code. 
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AGGREGATE MANIPOIATION PACKAGE 

The following routines constitute the Aggregate Manipulation Package 
(AMP). They support the CALL/360-OS built-in functions PROD, SUM, 
and POLY. Descriptions of these routines are given on succeeding pages 
of this manual, in the order indicated. 

Note that the four general groupings appear in alphabetic order 
according to the mnemonics of the routines. Within each grouping, 
discussion of short floating-point precedes discussion of long floating- 
point, and discussion of real precedes discussion of complex. 

Indexing: 

Interleaved Array Indexer (IHEJXI) 

PROD: 



PROD-Interleaved 
PROD- Inter 1 ea ved 
PROD-Interleaved 
PROD- In t er 1 ea ved 
PROD-Interleaved 
PROD-Interleaved 



Real Fixed Array (IHEPDF) 
Real Short Float Array (IHEPDS) 
Real Long Float Array (IHEPDL) 
Complex Fixed Array (IHEPDX) 
Complex Short Float Array (IHEPDW) 
Complex Long Float Array (IHEPDZ) 



SUM: 



SUM-Inter leaved Real Fixed Array (IHESMF) 
SUM-Interleaved Real/Complex Short Float Array (IHESMG) 
SUM-Inter leaved Real/Complex Long Float Array (IHESMH) 
SUM-Interleaved Complex Fixed Array (IHESMX) 



POLY: 



POLY 


(A,X) 


(A 


and 


X 


POLY 


(A,X) 


(A 


and 


X 


POLY 


(A,X) 


(A 


and 


X 


POLY 


(A,X) 


(A 


and 


X 


POLY 


(A,X) 


(A 


and 


X 


POLY 


(A,X) 


(A 


and 


X 



Real Fixed) (IHEYGF) 
Real Short Float) (IHEYGS) 
Real Long Float) (IHEYGL) 
Complex Fixed) (IHEYGX) 
Complex Short Float) (IHEYGW) 
Complex Long Float) (IHEYGZ) 



SPEED 

The average execution times given in this subsection are based on 
information in IBM System/360 Instruction Timing Information (A22-6825). 

EFFECT OF HEXADECIMAL TRUNCATION 

Allowance must be made for a certain amount of truncation error in 
System/360 handling of hexadecimal numbers. The accuracy of values 
returned by the built-in fimctions SUM, PROD, and POLY is governed 
accordingly. (See "Hexadecimal Truncation Errors'" under "Mathematical 
Function Package.") 

ARGUMENTS 

Any restrictions on arguments are noted under two headings: 

Range : This states any range of arguments for which a module is valid. 
Arguments outside the given ranges are assumed to have been excluded 
before the module is called. 
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Error and Except-ional Condi t:ions ; These cover conditions which may 
result from the use of a routine; they are listed in four categories: 

P - Programmed conditions in the module concerned • Programmed tests 
are made where not too costly and, if an invalid argument is found, 
a branch is taken to the entry point IHEERRB of the Error Routine, 
(See "Handling of Interrupts Package" in Volume II. Error Routine 
is part of EXEP, which is a subpackage of HIP,) An appropriate 
message is printed and the ERROR condition is raised, 

I - Interrupt conditions in the module concerned. For those routines 
where FIXEDOVERFLOW is detected by programmed tests, or where 
hardware interruptions may occur, the OVERFLOW, UNDERFLOW, and (when 
the conversion package is called) ERROR conditions pass to IHEERR 
and are treated in the normal way. The machine is assumed to be 
enabled for all interruptions except significance, which is masked. 

O - Programmed conditions in modules called by the module concerned. 
These occur when invalid arguments are detected in the module called. 

H - As I, but the interrupt conditions occur in the modules called 
by the module concerned. 
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TITLE: INTERLEAVED ARRAY INDEXER (IHEJXI) 

Program Definition 

Purpose and Usage 

Interleaved Array Indexer is used to find the next element of an array 
and to return its byte address in register P7. 

Entry point IHEJXIY is used to initialize the routine for byte addresses 
and to provide the address of the first element in the array. Entry 
point IHEJXI A is used thereafter to obtain the addresses of subsequent 
elements of the array? one address is returned for each entry into 
the routine. 

Description 

Method ; 

Arrays are stored in row-major order. Let Li be the lower bound and 
Ui the upper bound of the ith dimension, and n the number of dimensions. 
Starting with the element A(L1,L2,,. . . ,Ln), the routine varies the 
subscripts through their ranges to A(U1,U2, . . . ,Un),^ changing the nth 
subscript most rapidly; in this way the elements are referenced in 
the order in vi^ich they are stored. 

In the initializer part of the routine, it is determined how many of 
the last dimensions in the array are stored contiguously. This is 
done by comparing Mj for equality with U(j+1) - L(j+1) ♦ 1 for j = 
n-1, n-2,...„l. The contiguous storage of elements is ended at the 
first point where this equality does not hold. This dimension is saved 
in the storage location called S. The count of the number of contiguous 
elements is calculated during this process and saved. An array INDEX 
is set to the lower bound values for all dimensions. 

Initially the base address is equal to the address of A(L1,L2, . . . ,Ln) . 
Each subsequent contiguous element address is generated from the 
previous one by adding the multiplier Mn from the array dope vector 
(ADV) and reducing the contiguous element count by 1. 

When the count reaches 0, the next subscript to the left is increased 
by 1, the next element address is recalculated using the values in 
INDEX, and the contiguous element count is restored. If increasing 
the INDEX value for a dimension causes it to pass the upper bound 
value, then the next dimension to the left is increased by 1, and so 
forth. In any case,, when increasing a dimension by 1, all dimensions 
to the right are reset to their lower botmds. 

Range ; 

< number of dimensions < 255 
Implementation ; 

• Module size: 1196 bytes 
Errors Detected 
None 
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Local Variables 

NDIMS Niimber of dimensions 

ELE Niimber of contiguous elements 

COUNT Number of contiguous elements left 

Q Current dimension processing 

S Dimension ruining continuity 

MN Last multiplier 

BNDPTR Pointer to bound slot for dimension Q 

R Leftmost dimension being worked on 

UBNDPTR Pointer to upper bound for dimension Q 

INDEX Current subscript value array 

Note: Storage for these local variables is kept in the routine itself; 
thus this routine is not reentrant. However, it can easily 
be made reentrant by moving its storage to the LCA area, since 
the storage needs to be preserved between calls to this routine. 

Program Interface 



Entry Points 
IHEJXIY 



IHEJXIA 



Entry to initialize for byte addressing, and to locate 
the first element of the array. 

P7 = address of parameter list 

parameter list: A(Source ADV) 

Entry to locate the next element of the array 

Linkage: No explicit argiaments. 
Implicit arguments: LCA 

Exit Conditions 

Normal. Returns to caller via the link register with: 

1. P7 = the byte address of the first or current element of the array, 

2. GO = 0, if the last address of the last element of the array 

provided was not the last element. 

3. GO = 4, if the last address of the last element of the array 

provided was the last element. 

Routines Called 



None 




Global Variables 


WJXIDVA 


Dope V< 


WJXILADD 


Last e: 


Comments 




Called by: 




IHEPDF 


IHESMF 


IHEPDX 


IHESMX 


IHEPDS 


IHKSMG 


IHEPDW 


IHESMH 


IHEPDL 




IHEPDZ 
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TITLE: PROD- INTERLEAVED REAL FIXED ARRAY (IHEPDF) 

Program Definition 

Purpose and Usage 

PROD- Inter leaved Real Fixed Array is used to equate a long or short 
floating-point real target to the product of all the elements of an 
interleaved array of fixed-point reals. 

Description 

Method: 

The elements of the array are used in row-raajor order to multiply the 
current product. For fixed-point arguments, each element is converted 
to floating-point by using the Total Conversion Package. The precision 
specified in the source DED determines the precision of the target. 

Implementation : 

• Module size: 144 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the following formulas. 
■Short" or 'long* refers to the floating-point result. 

Target 

short a ♦ Tl ♦ R*(e + T3) 

long c + Tl + R*(f + T3) 

Constants used in the formulas are: 

R = number of elements in the array 

Tl = sum of times required to execute IHEJXI using IHEJXIY 
and IHEJXIA 

T3 = time for the appropriate conversion using IHEDMA 



r 




1 




IBM System/ 360 Model 


Number 








1 


30 


1 


HO 


1 


50 


1 


65 


1 


75 




a 


1 


1075 


1 


365 


1 


141 


1 


37.8 


1 


25 




c 


1 


112a 


1 


382 


1 


147 


1 


39 


1 


25.9 




e 


1 


eas 


1 


201 


1 


69 


1 


18.2 


1 


12.1 




f 


1 


1363 


1 


371 


1 


81.5 


1 


20.2 


1 


13 



Errors Detected 

Error and Exceptional Conditions : 

I: OVERFLOW (300) 
UNDERFLOW (340) 

Local Variables 

None 
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Program Interface 

Entry Points 

IHEPDFO 

P7 = A (Parameter List) where Parameter List: 
A(ADV) 

A (Array DED) 
A (Target) 

Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

IHEDMA Arithmetic Conversion Director 
IHEJXI Interleaved Array Indexer 

Global Variables 

None 

Comments 

Called by compiled code. 
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TITLE: PROD- INTERLEAVED REAL SHORT FLOAT ARRAY (IHEPDS) 

Program DefinitJ.on 

Purpose and Usage 

PROD-Inter leaved Real Short Float Array is used to ecjuate a short 
floating-point real target to the prodxict of all the elements of an 
interleaved array of short floating-point reals. 

Description 

Method ; 

The elements of the array are used in row-major order to multiply the 
current product. 

Implementation ; 

• Module size: 88 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the following formula. 

a + Tl + R*b 

Constants used in the formula are: 

R = number of elements in the array 

Tl = sum of times required to execute IHEJXI using IHEJXIY 
and IHEJXIA 



r 




1 




IBM System/ 360 Model Number 








1 


30 


1 «*0 1 50 1 65 1 


75 




a 
b 


1 
1 


ai8 

U92 


1 178 1 80 1 23.7 1 
1 11*3 1 43.3 1 11.4 1 


17.5 
7.2 



Errors Detected 

Error and Exceptional Conditions ; 

I; OVERFLOW (300) 
UNDERFLOW (340) 

Local Variables 

None 

Program Interface 

Entry Points 

IHEPDSO 

P7 = A (Parameter List) where Parameter List; 
A(ADV) 
A (Target) 
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Exit Conditions 

Normal. Returns to caller via the link register • 

Routines Called 

IHEJXI Interleaved Array Indexer 
Global Variables 
None 

Comments 
Called by compiled code. 
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TITLE: PROD- INTERLEAVED REAL LONG FLOAT ARRAY (IHEPDL) 

Program Definiiiion 

Purpose and Usage 

PROD-Inter leaved Real Long Float Array is used to equate a long 
floating-point real target to the product of all the elements of an 
interleaved array of long floating-point reals. 

Description 

Method: 

The elements of the array are used in row-major order to multiply the 
current product. 

Implementation : 

• Module size: 88 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the following formula. 

a ♦ Tl ♦ R*b 

Constants used in the formula are: 

R = number of elements in the array 

Tl = sum of times required to execute IHEJXI using IHEJXIY 
and IHEJXIA 



r 




1 




IBM System/360 Model Number 








1 


30 


1 40 1 50 1 65 1 


75 


1 


a 
b 


1 
! 


-290 
126it 


1 8.8 1 67.5 1 20.7 | 
I 331 1 63.8 1 m.8 1 


15.5 
9.2 



Errors Detected 

Error and Exceptional Conditions : 

I: OVERFLOW (300) 
UNDERFLOW (3«*0) 

Local Variables 

None 

Program Interface 

Entry Points 

IHEPDLO 

P7 = A (Parameter List) where Parameter List: 
A(ADV) 
A (Target) 
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Exit Conditions 

Normal, Returns to caller via the link register. 

Routines Called 

IHEJXI Interleaved Array Indexer 
Global Variables 
None 

Comments 
Called by compiled code. 
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TITLE; PROD- INTERLEAVED COMPLEX FIXED ARRAY (IHEPDX) 

Program Definite- on 

Purpose and Usage 

PROD-Inter leaved Complex Fixed Array is used to equate a long or short 
floating-point complex target to the product of all the elements of 
an interleaved array of fixed- point complexes. 

Description 

Method ; 

The elements of the array are used in row-major order to multiply the 
current product. For fixed-point arguments, each element is converted 
to floating-point by using the Total Conversion Package. The precision 
specified in the source DED determines the precision of the target. 

Implementation ; 

• Module size; 272 bytes 

• Execution times; 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the following formula. 



Source 



Target 
short 



binary short a ■•- Tl + R* (e + 2*T3) 

Constants used in the formula are; 

R = number of elements in the array 

Tl = sum of times required to execute IHEJXI using IHEJXIY 
and IHEJXIA 

T3 = time for the appropriate conversion using IHEDMA 



r 




1 




IBM Systero/360 Model 


Number 








1 


30 


1 


HO 


1 


50 


1 


65 


1 


75 




a 


1 


1382 


1 


H97 


1 


208 


1 


56.8 


1 


ai.2 




c 


1 


137U 


1 


501 


1 


205 


1 


56.7 


1 


JJ0.8 




e 


1 


201*7 


1 


574 


1 


182 


1 


fXH.H 


1 


26.3 




f 


1 


5097 


1 


1306 


1 


250 


1 


56.9 


1 


33.2 



Errors Detected 

Error and Exceptional Conditions; 

I; OVERFLOW (300) 
UNDERFLOW (340) 

Local Variables 

None 
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Program Interface 

Entry Points 

IHEPDXO 

P7 = ACParameter List) where Parameter List: 
A(ADV) 

A (Array DED) 
A (Target) 

Exit Conditions 

Normal* Returns to caller via the link register. 

Routines Called 

IHEDMA Arithmetic Conversion Director 
IHEJXI Interleaved Array Indexer 

Global Variables 

None 

Comments 

Called by compiled code. 
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TITLE: PROD-INTERLEAVED COMPLEX SHORT FLOAT ARRAY (IHEPDW) 

Program Definition 

Purpose and Usage 

PROD-Interleaved Complex Short Float Array is used to equate a short 
floating-point complex target to the product of all the elements of 
an interleaved array of short floating-point complexes. 

Description 

Method ; 

The elements of the array are used in row-major order to multiply the 
current product. 

Implementation ; 

• Module size: 120 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the following formula. 

a + Tl ♦ R*b 

Constants used in the formula are: 

R = number of elements in the array 

Tl = sum of times required to execute IHEJXI using IHEJXIY 
and IHEJXI A 



r 


1 




IBM Systein/360 Model Number 






1 


30 


1 40 


1 50 


1 65 1 


75 


1 ^ 
1 b 


1 
1 


-81U 
169it 


1 -96 
1 452 


\ 12.1 
1 132 


1 8.9 1 
1 32 1 


12.5 
17.7 



Errors Detected 

Error and Exceptional Conditions ; 

I: OVERFLOW (300) 
UNDERFLOW (340) 

Local Variables 

None 

Program Interface 

Entry Points 

IHEPDWO 

P7 = A (Parameter List) where Parameter List: 
A(ADV) 
A (Target) 
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Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

IHEJXI Interleaved Array Indexer 
Global Variables 
None 

Comments 
Called by compiled code. 
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TITLE; PROD- INTERLEAVED COMPLEX LONG FLOAT ARRAY (IHEPDZ) 

Program DefinitdLon 

Purpose and Usage 

PROD- Interleaved Complex Long Float Array is used to equate a long 
floating-point complex target to the product of all the elements of 
an interleaved array of long floating-point complexes. 

Description 

Method ; 

The elements of the array are used in row-major order to multiply the 
current product. 

Implementation ; 

• Module size; 120 bytes 

• Execution times; 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the following formula. 

a + Tl + R*b 

Constants used in the formula are: 

R = number of elements in the array 

Tl = sum of times required to execute IHEJXI using IHEJXIY 
and IHEJXIA 



r 




1 


IBM System/ 


360 Model 


Number 






1 


30 ! 


40 


1 


50 1 


65 1 75 




a 
b 


1 

1 


-38IH 1 
4830 I 


-852 
1214 


1 
1 


-61.5 1 
211 1 


-4.7 1 4.2 
46.1 1 25.7 



L 

Errors Detected 

Error and Exceptional Conditions ; 

I; OVERFLOW (300) 
UNDERFLOW (3t»0) 

Local Variables 

None 

Program Interface 

Entry Points 

IHEPDZ 

P7 = A (Parameter List) where Parameter List; 
A(ADV) 
A (Target) 
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Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

IHEJXI Interleaved Array Indexer 
Global Variables 
None 
Comments 
Called by compiled code. 
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TITLE: SUM- INTERLEAVED REAL FIXED ARRAY (IHESMF) 

Program Definition 

Purpose and Usage 

SUM-Interleaved Real Fixed Array is used to equate a long or short 
floating-point real target to the stim of all the elements of an 
interleaved array of fixed- point reals. 

Description 

Method : 

The elements of the array are added to the current sum in row-major 
order. For fixed-point arguments, each element is converted to 
floating-point by using routines of the Total Conversion Package. 
The precision specified in the source DED determines the precision 
of the target. 

Implementation : 

• Module size: 136 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the following formulas. 
•Short" or 'long' refers to the floating-point result. 

Target 

short a -t- Tl ♦ R*(e ♦ T3) 

long c ♦ Tl ♦ R*(f + T3) 

Constants used in the formulas are: 

R = number of elements in the array 

Tl = sum of times required to execute IHEJXI using IHEJXIY 
and IHEJXIA 

T3 = time for the appropriate conversion using IHEDMA 



r 




1 




IBM System/ 360 Model 


Number 








1 


30 


1 


40 


1 


50 




1 


65 1 


75 




a 


1 


107U 


1 


363 


1 


141 




1 


37.3 1 


24.6 




c 


1 


1123 


1 


381 


1 


147 




1 


38.5 1 


25.6 




e 


1 


400 


1 


139 


1 


54, 


4 


1 


16.2 1 


10.9 




f 


1 


U18 


1 


139 


1 


53, 


2 


1 


15.1 1 


9.8 



Errors Detected 

Error and Exceptional Conditions : 

I: OVERFLOW (300) 
UNDERFLOW (340) 
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Local Variables 

None 

Program Int:erface 

Entry Points 

IHESMFO 

P7 = A (Parameter List) where Parameter List: 
A(ADV) 

A (Array DED) 
A (Target) 

Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

IHEDMA Arithmetic Conversion Director 
IHEJXI Interleaved Array Indexer 

Global Variables 

None 

Comments 

Called by compiled code. 
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TITLE: SUM- INTERLEAVED REAL/COMPLEX SHORT FLOAT ARRAY (IHESMG) 

Program Definii:ion 

Purpose and Usage 

SUM-Interleaved Real/Complex Short Float Array is used to equate a 
short floating-point real or complex target to the sum of all the 
elements of an interleaved array of short floating-point reals or 
complexes , respectively. 

Description 

Method ; 

The elements of the array are added to the current sum in row-major 
order. For a complex argument, the summations of the real and imaginary 
parts are developed concurrently. 

Implementation ; 

• Module size: 128 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the following formulas. 

Real a + Tl ♦ R*b 

Complex c + Tl ♦ R*d 

Constants used in the formulas are: 

R = number of elements in the array 

Tl = sum of times required to execute IHEJXI using the IHEJXIY 
and IHEJXIA 



r 




1 




IBM System/ 360 Model 


Number 


1 
1 






1 


30 


1 


40 1 


50 1 


65 


1 


1 

75 1 




a 
b 
c 
d 


1 
1 
1 
1 


935 

167 

1129 

334 


1 
1 
1 
1 


314 1 
43.9 1 

372 1 
87.8 1 


121 1 
15.2 1 

142 1 
30.4 1 


32 

4.2 
36.8 

8.3 


1 
1 
1 
1 


21.4 1 

2.1 1 
24.6 1 

4.2 1 



Errors Detected 

Error and Exceptional Conditions : 

I: OVERFLOW (300) 
UNDERFLOW (340) 

Local Variables 

None 
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Program Interface 

Entry Points 

IHESMGR Entry for real arrays. 

P7 = A(Parameter List) where Parameter List: 
A(ADV) 
A (Target) 

IHESMGC Entry for complex arrays. 
Linkage as for IHESMGR 
Exit Conditions 

Normal. Returns to caller via the link register. 
Routines Called 

IHEJXI Interleaved Array Indexer 
Global Variables 
None 
Comments 
Called by compiled code* 
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TITLE: SUM-INTERLEAVED REAL/COMPLEX LONG FLOAT ARRAY (IHESMH) 

Proqrain Definition 

Purpose and Usage 

SUM-Interleaved Real/Coroplex Long Float Array is used to equate a long 
floating-point real or complex target to the sum of all the elements 
of an interleaved array of long floating-point reals or complexes, 
respectively • 

Description 

Method ; 

The elements of the array are added to the current sum in row-major 
order. For a complex argument, the summations of the real and imaginary 
parts are developed concurrently. 

Implementation ; 

• Module size: 128 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the following formulas. 

Real a + Tl + R*b 

Complex c ♦ Tl ♦ R*d 

Constants used in these formulas are: 

R = number of elements in the array 

Tl = sum of times required to execute IHEJXI using IHEJXIY 
and IHEJXIA 



IBM System/360 Model Number 



30 



I 



i\0 



I 



50 



I 



65 



I 



75 



a 
b 
c 
d 



I 



887 
366 
993 

5m 



I 



298 I 111 I 

116 I U2.t* I 

337 I 125 I 

157 I 55.9 I 



30.9 I 20.4 

12.3 I 7.9 

3C».U I 23.5 

16 I 9.8 



Errors Detected 

Error and Exceptional Conditions : 

I: OVERFLOW (300) 
UNDERFLOW (340) 

Local Variables 

None 



107 



Program Interface 

Entry Points 

IHESMHR Entry for real arrays • 

P7 = A (Parameter List) where Parameter List: 
A(ADV) 
A (Target) 

IHEl^HC Entry for complex arrays. 
Linkage as for IHESMHR 
Exit Conditions 

Normal. Returns to caller via the link register. 
Routines Called 

IHEJXI Interleaved Array Indexer 
Global Variables 
None 

Comments 
Called by compiled code. 
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TITLE: SUM-INTERLEAVED COMPLEX FIXED ARRAY (IHESMX) 

Program Definition 

Purpose and Usage 

SUM-Interleaved Complex Fixed Array is used to equate a long or short 
floating-point complex target to the sum of all the elements of an 
interleaved array of fixed- point complexes. 

Description 

Method: 

The elements of the array are added to the current sum in row-major 
order. For fixed-point arguments, each element is converted to 
floating-point by using routines of the Total Conversion Package. 
For a complex argument, the summations of the real and imaginaiy parts 
are developed concurrently. The precision specified in the source 
DED determines the precision of the target. 

Implementation : 

• Module size: 224 bytes 

• Execution times: 

Approximate execution times in microseconds for the System/360 
models given below are obtained from the following formula- 



Source 



Target 
short 



binary short a + Tl + R* (e + 2*T3) 

Constants used in the formula are: 

R = ntimber of elements in the array 

Tl = sum of times required to execute IHEJXI using IHEJXIY 
and IHEJXIA 

T3 = time for the appropriate conversion using IHEDMA 



r 




1 




IBM 


System/ 360 Model 


Number 








1 


30 


1 


40 


1 


50 


1 


65 1 


75 




a 


1 


1327 


1 


511 


1 


229 


1 


63.8 1 


a7.8 




c 


1 


1319 


1 


515 


1 


227 


1 


63.7 1 


m.n 




e 


1 


712 


1 


232 


1 


87 


1 


25.1 1 


16.4 




f 


1 


770 


! 


240 


1 


89 


1 


2U 1 


15.2 



Errors Detected 

Error and Exceptional Conditions : 

I: OVERFLOW (300) 
UNDERFLOW (340) 

Local Variables 

None 
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Program Interface 

Entry Points 

IHESMXO 

P7 = ACParameter List) where Parameter List: 
A(ADV) 

ACDED for Array) 
A (Target) 

Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

IHEDMA Arithmetic Conversion Director 
IHEJXI Interleaved Array Indexer 

Global Variables 

None 

Comments 

Called by compiled code. 
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TITLE: POLY (A, X) (A AND X REAL FIXED) (IHEYGF) 

Program Definition 

Purpose and Usage 

POLY (A,X) for real fixed A and X has two options: 

1. Vector X: 

Let the argiiments be arrays declared as A(m:n) and X(p:q). 
Then the function computed is: 

n - 111 j - 1 
A(in) + \^ A(m +j)'«^ J^ X(p + i) 
i= 1 i==0 

unless n = m^ when result is A(ro) . 
If 

q-p<n-m-l, 
then t for 

p + i > q, X(p+i) = X(q), 

2. Scalar X: 

This may be interpreted as a special case of vector X, that 
is, a vector with one element, X(l), which is equal to X, Then 
the function computed is: 

n - m 

y ^ A(m +j)'J=X=i i-| 
1 = 

A floating-point result is obtained in both cases. 

Description 

Method : 

1. Vector X, (q - p > n - m - 1): 

POLY (A,X> is evaluated by nested multiplication and addition, 

that is, (... (A(n)*X(k)+A(n-l))*X(k-l)+A(n-2))*...+A(m+l))*X(p)+A(m) 

where k=p+n-m-l. 

2. Vector X, (q - p < n - m - 1) : 

In the expression above, the terms in X with subscript ranging 
from k down to q are all made equal to X(q). The evaluation 
is treated as for scalar X until sufficient terms in X have 
been made equal to X(q), when the computation continues as in 1. 

3. Scalar X: 

Terms in X with subscript ranging from k to p are equal to X. 

For fixed-point arguments^ each element is converted to floating- 
point, by using the Total Conversion Package. 
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4, The target: precision is determined according to the highest 
precision specified in the source DED's. 

Implementation ; 

• Module size: 432 bytes 

• Execution times: 

Let the arguments be declared as A(m:n) and X(p:q)y or X, and T 
be the time for one conversion using the Arithmetic Conversion 
Director (IHEDMA) . Then the approximate execution times in 
microseconds for the System/360 models given are obtained from 
the appropriate formula. 'Short' or 'long' refers to the floating- 
point result. 

Scalar X: 

short a ♦ 2*T -•• (n-m) * (b+T) 
long c -I- 2*T + (n-m) ♦ (d+T) 

Vector X^ (q - p > n - m - 1): 



short 
long 



e ♦ T + 
g + T + 



(n-m) 
(n-m) 



(f ♦ 2*T) 
(h + 2*T) 



Vector X, (q - p < n - m - 1) : 



short i ♦ 2*T + (n-m) ♦ (b+T) + (q-p+1) ♦ (j+T) 
long k ♦ 2*T + (n-ro) * (d+T) + (q-p+1) * (1+T) 



IBM System/360 Model Number 



30 



I 



40 



50 



I 



65 



75 



a 
b 
c 
d 
e 
f 

g 

h 
i 

J 

k 

1 



2297 1 


834 


1 338 


1 99 


67.1 


904 


387 


1 100 


1 26 


17.1 


2408 


870 


351 


1 102 


69.2 


1706 


484 


124 


1 30.4 


20.2 


2629 


910 


370 


1 97.1 


64.6 


1480 


459 


155 


1 41.2 


1 26.3 


2740 1 


946 


383 


99.4 


66.7 


3044 


844 


1 200 


1 50 


1 32.6 


3370 


1197 


491 


140 


94 


258 


1 91 


35 1 


9.5 1 


6 


3581 1 


1233 


1 504 


143 


96.1 


280 


101 


1 39 


1 10.7 


7.1 



Errors Detected 

Error and Exceptional Conditions : 

I: OVERFLOW (300) 
UNDERFLOW (340) 

Local Variables 

None 
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Program Interface 
Entry Points 
IHEYGFV 
P7 



Entry for vector A and vector X. 



IHEYGFS 



= A (Parameter List) where Parameter List: 
A (A ADV) 
A(A DED) 
A(X ADV) 
A(X DED) 
A (Target) 

Entry for vector A and scalar X. 



P7 


A(PLIST) 


PLIST 


A (A ADV) 




A(A DED) 




A(X) 




A(X DED) 




A (Target) 



Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

IHEDMA Arithmetic Conversion Director 
Global Variables 
None 
Comments 
Called by compiled code. 
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TITLES POLY (A,X) (A AND X REAL SHORT FLOAT) (IHEYGS) 

Program Definition 

Purpose and Usage 

POLY (A^X) for real short float A and X has two options: 

1. Vector X: 

Let the arguments be arrays declared as A(m:n) and X(p:q). 
Then the function computed is: 





n - m j - 1 
A(ni) + V^ A(m +j) * J[J[ X(p + i) 

j = 1 i = 




unless n = m, when result is A(m) 




If 




q-p<n-ro-l. 




then^ for 




p + i > q, X(p+i) = X(q). 


2. 


Scalar X: 



This may be interpreted as a special case of vector X^ that 
is, a vector with one element, X(l), which is equal to X. Then 
the function computed is: 

n - m 

y A(m + j)*X**j 
j = 

A floating-point result is obtained in both cases • 

Description 

Method: 

1. Vector X, (q - p > n - m - 1): 

POLY (A,X) is evaluated by nested multiplication and addition, 

that is, (... (A(n)*X(k)+A(n-l))*X(k-l)+A(n-2))*..,+A(m+l))*X(p)+A(m) 

where k=p+n-m-l. 

2. Vector X, (q - p < n - m - 1) : 

In the expression above, the terms in X with subscript ranging 
from k down to q are all made equal to X(q). The evaluation 
is treated as for scalar X until sufficient terms in X have 
been made equal to X(q), when the computation continues as in 1. 

3. Scalar X: 

Terms in X with subscript ranging from k to p are equal to X, 

4. The target is short precision, 



114 



Implementation ; 

• Module size: 2^0 bytes 

• Execution times: 

Let the arguments be declared as A(m:n) and X(p:q), or X, and T 
be the time for one conversion using the Arithmetic Conversion 
Director (IHEDMA). Then the approximate execution times in 
microseconds for the System/360 models given are obtained from 
the appropriate formula. 



Scalar X: 








a ♦ (n - m) *b 








Vector X: 








(q-p>n-m- 


1): 


c + (n - 


- m)*d 


(q-p<n-m- 


1): 


e + (n - 


- m)*b + 



(q - p ♦ l)*f 



a 
b 
c 
d 
e 
f 



I 



IBM System/360 Model Ntimber 



30 



I 



i\0 



I 



50 



65 



I 



75 



1232 


«^30 


182 


49.5 


1 33.3 


461 


121 


37.6 


9.8 


1 5.3 


1871 


623 


259 


69.6 


1 45.4 


U90 


128 


40.9 


10.5 


1 5.7 


21U0 


733 


304 


83.6 


1 54.7 


29 


7.5 


3.3 


0.7 


0.4 



L J 

Errors Detected 

Error and Exceptional Conditions ; 

I: OVERFLOW (300) 
UNDERFLOW (340) 



Local Variables 


None 




Proqram Interface 


Entry Points 




IHEYGSV 


Entry for vector A and vector X 


P7 
PLIST 


A (PLIST) 
A (A ADV) 
A(X ADV) 
A (Target) 


IHEYGSS 


Entry for vector A and scalar X 


P7 


= A(Parameter List) where Parameter List: 
A (A ADV) 
A(X) 
A (Target) 
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Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

None 

Global Variables 

None 

Comments 

Called by compiled code. 
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TITLE: POLY (A.,X) (A AND X REAL LONG FLOAT) (IHEYGL) 

Program Definition 

Purpose and Usage 

POLY (A,X) for real long float A and X has two options: 

1. Vector X: 

Let the arguments be arrays declared as A(m:n) and X(p:q). 
Then the function computed is: 

n-rn j - 1 

A(m)+ y. A(m+j)*J][ X(p + i) 

j= 1 i = 

unless n = m^ when result is A(m) . 
If 

q-p<n-m-l,^ 
then^ for 

p ♦ i > q^ X(p+i) = X(q). 

2. Scalar X: 

This may be interpreted as a special case of vector X, that 
is,, a vector with one element, X(l), which is equal to X. Then 
the function computed is: 

n - m 

/ A(m + j)*X**j 

.1 = 

A floating-point result is obtained in both cases. 

Description 

Method: 

!• Vector X, Cq - p > n - m - 1): 

POLY CA,X) is evaluated by nested multiplication and addition, 

that is, (... CA(n)*X(k)+A(n-l))*X(k-l)+A(n-2))*...+A(m+l))*X(p)+ACm] 

where k=p+n-m-l. 

2. Vector X„ (q - p < n - m - 1) : 

In the expression above, the terms in X with subscript ranging 
from k down to q are all made equal to XCq), The evaluation 
is treated as for scalar X until sufficient terms in X have 
been made equal to X(q), when the computation continues as in 1. 

3. Scalar X: 

Terms in X with subscript ranging from k to p are equal to X. 

4. The target is long precision. 
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Iinplemeni:at:ion : 

• Module size: 240 bytes 

• Execution times: 

Let the arguments be declared as A(m:n) and X{p:q), or X, and T 
be the time for one conversion using the Arithmetic Conversion 
Director (IHEDMA). Then the approximate execution times in 
microseconds for the System/360 models given are obtained from 
the appropriate formula. 



Scalar X: 




a + (n - m)*b 




Vector X: 




(q ~ p > n - m - 
(q-p<n-m- 


- 1) 

- 1) 



c + (n - m)*d 

e + (n - m)*b + (q - p ♦ 1) *f 



30 



IBM System/360 Model Number 



I 



40 



50 



65 



I 



75 



a 
b 
c 

d 
e 

f 



I 



1320 


1 «t51 


1 189 


49.6 


1 33.3 


12m 


308 


1 56.9 


13.1 


7.3 


1959 


644 


1 266 


69.7 


45.4 


1270 


1 316 


1 60.2 


13.7 


7.7 


2228 


1 755 


1 311 


83.6 


54.7 


29 


7.5 


3.3 


0.7 1 


0.4 



L J 

Errors Detected 

Error and Exceptional Conditions ; 

I: OVERFLOW (300) 
UNDERFLOW (3^0) 



Local Variabl 


es 


None 




Proqram Interface 


Entry Points 




IHEYGLV 


Entry for vector A and vector X, 


P7 


A (PLIST) 


PLIST 


A (A ADV) 




A(X ADV) 




A (Target) 


IHEYGLS 


Entry for vector A and scalar X. 


P7 = 


A (Parameter List) where Parameter List: 




A (A ADV) 




A(X) 




A (Target) 
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Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

None 

Global Variables 

None 

Comments 

Called by compiled code. 
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TITLE: POLY (A^X) (A AND X COMPLEX FIXED) (IHEYGX) 

Program Definit:ion 

Purpose and Usage 

POLY (A^X), for complex fixed A and X has two options: 

1. Vector X: 

Let the arguments be arrays declared as A(m:n) and X(p:q). 
Then the function computed is: 



n - m j - 1 

A(m) + N A(m + j)*TT 

J = 1 i = 



unless n = m^ when result is A(m). 

If 

q-p<n-m-l^ 

then , for 

p + i > q, X(p+i) = X(q). 

2, Scalar X: 

This may be interpreted as a special case of vector X, that 
is^ a vector with one element,, X(l), which is equal to X. Then 
the function computed is: 

n - m 

y A(m + j)*X**j 

j = 

A floating-point result is obtained in both cases. 

Description 

Method : 

!• Vector X„ (q - p > n - m - 1) : 

POLY (AfX) is evaluated by nested multiplication and addition, 

that is, (... (A(n)*X(k)+A(n-l>)*X(k-l)+A(n-2))*...+A(m+l))*X(p)+A(m) 

where k=p+n-m-l. 

2. Vector X, Cq - p < n - m - 1) : 

In the expression above, the terms in X with subscript ranging 
from k down to q are all made equal to X(q). The evaluation 
is treated as for scalar X until sufficient terms in X have 
been made equal to X(q), when the computation continues as in 1. 

3. Scalar X: 

Terms in X with subscript ranging from k to p are equal to X, 

For fixed- point arguments, each element is converted to floating- 
point, by using the Total Conversion Package. 
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4. 



The target precision is determined by the highest precision 
indicated in the array/scalar DED's. 



Implementation ; 

• Module size: 688 bytes 

• Execution times: 

Let the arguments be declared as A(m:n) and X(p:q), or X, and T 
be the time for one conversion using the Arithmetic Conversion 
Director (IHEDMA) . Then the approximate execution times in micro- 
seconds for the System/360 models given are obtained from the ap- 
propriate formula. 'Short' or 'long* refers to the floating-point 
result. 



Scalar X: 

short 
long 



a + 4*T + (n-m)*(b+2*T) 
c ♦ i^*T + (n-m)*(d+2*T) 



Vector X, (q - p > n - m - 1) 



short 
long 



e + 2*T + (n-m)*(f+4*T) 
g + 2*T + (n-m)*(h+4*T) 



Vector X, (q-p<n-m-l): 



short 
long 



i + u*T + (n-m)*(b+2*T)+(q-p+l)*(j+2*T) 
k + 4*T + (n-m)*(d+2*T)+(q-p+l)*(l+2*T) 



r 




1 IBM System/360 Model 


Number 








1 30 


1 40 


1 50 


1 


65 


1 75 




a 


1 32U5 


1 1174 


1 478 




142 


1 95.4 




b 


1 23U5 


1 664 


1 220 




58.3 


1 33.5 




c 


1 3471 


1 1221 


496 




145 


1 97.6 




d 


1 5519 


1 1433 


1 301 




73.2 


1 42.6 




e 


1 3187 


1 1114 


447 




127 


1 83.5 




f 


1 4368 


1226 


399 




103 


1 58.1 




g 


1 3533 


1 1161 


465 




129 


1 85.7 




h 


1 10636 


2746 


556 




133 


1 76.4 




i 


1 4087 


1459 


592 




171 


1 114 




J 


1 545 


1 187 1 


72. 


3 1 


21.1 1 


12.7 




k 


1 4243 


1506 1 


610 




176 


1 116 




1 


1 567 


1 196 1 


76. 


3 1 


22.3 1 


13.6 



Errors Detected 

Error and Exceptional Conditions : 

I: OVERFLOW (300) 
UNDERFLOW (340) 

Local Variables 

None 
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Program Interface 
Entry Points 
IHEYGXV 
P7 



Entry for vector A and vector X, 



IHEYGXS 



= ACParcuneter List) where Parameter List: 
A (A ADV) 
A (A DED) 
A(X ADV) 
A(X DED) 
A (Target) 

Entry for vector A and scalar X. 



P7 


A(PLIST) 


PLIST 


A (A ADV) 




AC A DED) 




ACX) 




ACX DED) 




A (Target) 



Exit Conditions 

Normal, Returns to caller via the link register. 

Routines Called 

IHEDMA Arithmetic Conversion Director 
Global Variables 
None 
Comments 
Called by compiled code. 
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TITLE: POLY (A^X) (A AND X COMPLEX SHORT FLOAT) (IHEYGW) 

Program Definition 

Purpose and Usage 

POLY (A,X) for complex short float A and X has two options: 

1. Vector X: 

Let the arguments be arrays declared as A(m:n) and X(p:q). 
Then the function computed is: 

n - m j - 1 
A(m) + y^ A(m + j) * JJ X(p + i) 

j=l i=0 

unless n = m^ when result is A(m) • 
If 

q-p<n-m-l, 
then, for 

p + i > q, X(p+i) = X(q). 

2. Scalar X: 

This may be interpreted as a special case of vector X,, that 
iSt a vector with one element, X(l), which is equal to X, Then 
the function computed is: 

n - m 

y A(m + j)*X**j 

j = 

A floating-point result is obtained in both cases. 
Description 
Method : 

1. Vector X,, (q - p > n - m - 1): 

POLY (A,X) is evaluated by nested multiplication and addition, 

that is, (•.. (A(n)*X(k)+A(n-l))*X(k-l)+A(n-2))*---*A(m+l))*X(p)+A(m) 

where k = p-i-n-m-l« 

2. Vector X, (q - p < n - m - 1) : 

In the expression above, the terms in X with subscript ranging 
from k down to q are all made equal to XCq). The evaluation 
is treated as for scalar X until sufficient terms in X have 
been made equal to X(q), when the computation continues as in 1. 

3, Scalar X: 

Terms in X with subscript ranging from k to p are equal to X, 

4, The target is short precision. 
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Implementation ; 

• Module size: 280 bytes 

• Execution tiroes: 

Let the arguments be declared as A(m:n) and X(p:q)^ or X, and T 
be the time for one conversion using the Arithmetic Conversion 
Director (IHEDMA). Then the approximate execution times in micro- 
seconds for the System/360 models given are obtained from the 
appropriate formula. 

Scalar X: 



a ♦ (n ° ro)*b 






Vector X: 






(q-p>n-m- 


- 1): 


c + (n - m)*d 


(q-p<n-m- 


- 1): 


e + (n - m)*b 



♦ (q - p ♦ l)*f 



30 



IBM System/360 Model Number 



I 



t\0 



I 



50 



65 



I 



75 



a 
b 
c 
d 
e 
f 



I 



1396 


i»75 


198 


5l».9 1 


36.1 


1672 


1 i\25 


126 


30.5 1 


15 


2035 


667 


275 


75 1 


48.1 


1701 


1 a32 


129 


31.2 1 


15.4 


230i» 


775 


320 


88.9 1 


57.1* 


29 


1 7.5 


3.3 


0.7 1 


0.4 



L J 

Errors Detected 

Error and Exceptional Conditions : 

I: OVERFLOW (300) 
UNDERFLOW (340) 



Local Variables 


None 




Proqram Interface 


Entry Points 




IHEYGWV 


Entry for vector A and vector X 


P7 


A (PLIST) 


PLIST 


A (A ADV) 




A(X ADV) 




A (Target) 


IHEYGWS 


Entry for vector A and scalar X 


P7 


= A (Parameter List) where Parameter List: 




A (A ADV) 




A(X) 




A (Target) 
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Exit Conditions 

Normal. Returns to caller via the link register. 

Routines Called 

None 

Global Variables 

None 

Comments 

Called by compiled code* 
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TITLE: POLY (A,X) (A AND X COMPLEX LONG FLOAT) (IHEYGZ) 

Program Definite. on 

Purpose and Usage 

POLY (A,X), for complex long float A and X has two options: 

1, Vector X: 

Let the arguments be arrays declared as A(m:n) and X(p:q), 
Then the function computed is: 

n - m j - 1 

A(m) + y^ A(m + j)* J][ X(p + i) 

j= 1 i = 

unless n = m^ when result is A(m) . 
If 

q-p<n-m-lt 
then^ for 

p + i > q, X(p+i) = X(q). 

2. Scalar X: 

This may be interpreted as a special case of vector X, that 
is^ a vector with one element, X(l), which is equal to X. Then 
the function computed is: 

n - m 

y A(m + j)*X**j 

j = 

A floating-point result is obtained in both cases. 
Description 
Method ; 

1. Vector X,, (q - p > n - m - 1) : 

POLY (A,X) is evaluated by nested multiplication and addition, 

that is, (•..♦A(n)*X(k)+A(n-l))*X(k-l)+A(n-2))*...+A(m+l))*X(p)^A(m) 

where k=p+n-m-l. 

2. Vector X, (q - p < n - ro - 1) : 

In the expression above, the terms in X with subscript ranging 
from k down to q are all made equal to X(q), The evaluation 
is treated as for scalar X until sufficient terms in X have 
been made equal to X(q), when the computation continues as in 1. 

3. Scalar X: 

Terms in X with subscript ranging from k to p are equal to X. 

For fixed-point arguments, each element is converted to floating- 
point, by using the Total Conversion Package. 
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4. The target is long precision. 
Implementation ; 

• Module size: 280 bytes 

• Execution times: 

Let the arguments be declared as A(m:n) and X(p:q)t or X, and T 
be the time for one conversion using the Arithmetic Conversion 
Director (IHEDMA). Then the approximate execution times in 
microseconds for the Systero/360 models given are obtained from 
the appropriate formula. 



Scalar X: 




a + (n - m)*b 




Vector X: 




(q - p > n - m - 
(q-p<n~m- 


- 1) 

- 1) 



c + (n - m) *d 

e + (n - m)*b + (q - p + 1) ♦f 



a 
b 
c 
d 
e 
f 



I 



30 



IBM System/360 Model Number 



I 



UO 



I 



50 



65 



I 



75 



1572 


1 517 


1 211 


54.9 


36.1 


4824 


1183 


203 


44.2 


23 


2211 


1 710 


1 288 


75 


48.1 


4853 


1 1192 


207 


44.9 


23.4 


2480 


1 821 


1 333 


89 


57.4 


29 


7.5 


3.3 


0.7 


0.4 



L J 



Errors Detected 

Error and Exceptional Conditions : 

I: OVERFLOW (300) 
UNDERFLOW (340) 



Local Variabl 


es 


None 




Proqram Interface 


Entry Points 




IHEYGZV 


Entry for vector A and vector X. 


P7 
PLIST 


A (PLIST) 
A (A ADV) 
A(X ADV) 
A (Target) 


IHEYGZS 


Entry for vector A and scalar X, 


P7 


= A (Parameter List) where Parameter List: 
A (A ADV) 
A(X) 
A (Target) 
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Exit Conditions 

Normal, Returns to caller via the link register, 

Routines Called 

None 

Global Variables 

None 

Comroents 

Called by compiled code. 
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